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Abstract 

A  hybrid  method  of  moments  (MM)  based  numerical  model  for  the  electromagnetic  scattering 
from  large  finite  by  infinite  planar  slot  arrays  is  developed.  The  method  incorporates  the  novel 
concept  of  a  physical  basis  function  (PBF)  to  dramatically  reduce  the  number  of  required  unknowns. 
The  model  can  represent  a  finite  number  of  slot  columns  with  slots  oriented  along  the  infinite  axis, 
surrounded  by  an  arbitrary  number  of  coplanar  dielectric  slabs.  Each  slot  column  can  be  loaded  with  a 
complex  impedance,  allowing  one  to  tailor  the  edge  currents  to  provide  a  desired  echo  width  pattern. 

The  surface  equivalence  theorem  is  used  to  convert  the  original  slotted  ground  plane  geometry 
to  an  equivalent  unbroken  ground  plane  with  magnetic  surface  currents.  An  integral  equation  based 
on  these  magnetic  scattering  currents  is  solved  via  the  MM.  The  magnetic  currents  are  approximated 
by  a  set  of  basis  functions  composed  of  periodic  basis  functions  representing  the  edge  slot  columns 
and  a  single  PBF  representing  the  interior  slot  columns.  In  particular,  the  PBF  captures  the  behavior 
of  the  central  portion  of  the  array  where  the  perturbations  from  the  edges  have  become  negligible. 
Based  on  Floquet’s  theorem,  the  PBF  is  able  to  represent  an  arbitrarily  large  number  of  slot  columns 
with  just  one  unknown.  The  array  scanning  method  (ASM)  provides  the  contributions  from  the 
individual  edge  columns.  Finally,  a  newly  developed  one-sided  Poisson  sum  formulation  provides 
an  efficient  means  to  account  for  the  stratified  dielectric  media  via  a  spectral  domain  conversion. 

The  hybrid  method  is  validated  using  both  MM  reference  codes  and  measured  data.  The  results 
clearly  demonstrate  the  method’s  accuracy  as  well  as  its  ability  to  handle  array  problems  too  large  for 
traditional  MM  solutions. 
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7.  Introduction 

Most  of  the  fundamental  ideas  of  science  are  essentially  simple,  and  may,  as  a  rule,  be 

expressed  in  a  language  comprehensible  to  everyone. 

-  Albert  Einstein 

Periodic  screens  have  long  been  of  interest  to  the  scientific  and  engineering  communities.  Com¬ 
monly  called  frequency  surfaces  (FSS)  for  their  frequency  filtering  properties,  they  find  application 
over  a  large  range  of  the  electromagnetic  spectrum.  They  hold  particular  interest  to  the  Air  Force  as 
antenna  radomes  and  dichroic  antenna  feed  systems. 

FSS  typically  consist  of  regular  arrays  of  either  conducting  patch  elements  deposited  on  a 
dielectric  substrate  or  apertures  in  a  conducting  surface.  Figure  1  depicts  an  example  FSS  geometry. 
Arrays  of  this  type  exhibit  total  reflection  (patches)  or  transmission  (apertures)  near  the  resonance 
of  the  individual  elements.  In  studying  the  electromagnetic  scattering  from  these  structures,  most 
authors  use  an  integral  equation  formulation  which  connects  the  currents  induced  on  the  periodic 
elements  to  the  fields  scattered  from  the  FSS.  To  keep  the  problem  tractable,  they  approximate  the 
actual  curved  finite  array  with  a  flat  infinite  array  and  use  Floquet  analysis  to  reduce  the  infinite 
domain  problem  to  one  involving  a  representative  cell.  The  Method  of  Moments  (MM)  is  then  used 
to  convert  the  integral  equation  into  a  system  of  linear  equations  suitable  for  numerical  solution. 

Traditionally  the  actual  implementation  of  the  MM  in  infinite  array  analysis  has  taken  two 
different  tracks.  The  first,  pioneered  by  Munk  and  others  at  The  Ohio  State  University  [26,  28], 
approximates  the  induced  currents  as  filamentary  currents.  The  authors  use  a  plane  wave  expansion 
method  to  derive  very  efficient  expressions  for  the  radiated  fields.  The  Periodic  Moment  Method 
(PMM)  code  written  by  Henderson  [19]  is  an  excellent  example  of  the  use  of  this  technique. 

The  second  implementation,  developed  by  Mittra  and  others  at  the  University  of  Illinois  [23], 
favors  more  complex  basis  functions  to  better  model  the  actual  element  currents.  In  this  case  the 
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Figure  1.  Typical  FSS  Geometry 


authors  transform  the  integral  equation  to  the  spectral  domain  where  they  use  a  conjugate  gradient 
technique  to  solve  the  system.  Although  generally  more  accurate  for  a  given  number  of  basis 
functions,  codes  based  on  this  approach  tend  to  be  less  efficient. 

Regardless  of  the  solution  method,  all  of  the  infinite  array  techniques  fail  to  account  for  the 
edge  diffraction  or  curvature  effects  present  in  real  world  arrays.  A  knowledge  of  these  effects  is 
critical  to  advanced  radome  and  antenna  design.  Consequently,  recent  research  efforts  have  focused 
on  finite  array  analysis  methods. 

Due  to  their  finite  nature,  Floquet  theory  does  not  strictly  apply  to  these  arrays.  The  current 
on  each  of  the  array  elements  must  be  treated  as  an  independent  quantity.  Most  of  the  research  in 
this  area  falls  into  one  of  two  categories;  studying  the  effects  of  termination  using  relatively  small 
finite-by-infinite  or  finite-by-finite  arrays  or  finding  efficient  ways  to  compute  element  currents  in  an 
effort  to  handle  large  array  sizes.  In  both  cases,  the  general  approach  begins  with  the  development 
of  an  integral  equation  involving  unknown  element  currents.  Some  authors  choose  to  represent  the 
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integral  equation  kernel  in  the  spatial  domain  where  others  prefer  the  spectral  domain.  Next,  one 
approximates  the  currents  using  sub-domain  or  entire  domain  basis  functions.  With  an  appropriately 
defined  inner  product,  one  can  convert  the  integral  equation  to  a  matrix  equation  which  can  then  be 
solved  using  either  Gaussian  elimination  or  iterative  techniques. 

Examples  of  this  general  approach  include  Kasmer  and  Mittra’s  [20]  paper  based  on  the  spectral 
iteration  and  conjugate  gradient  iteration  techniques  applied  to  finite-by-finite  planar  patch  arrays. 
Grounds  and  Webb’s  [14]  attack  of  the  finite-by-infinite  problem  using  a  spectral  domain  Green’s 
function  formulation,  and  Skinner,  Munk,  and  Barre’s  [34,  2]  work  based  on  the  array  scanning 
method  (ASM)  applied  to  the  finite-by-infinite  problem  in  a  layered  dielectric  media.  All  of  these 
approaches  are  able  to  explicitly  model  the  edge  termination  effects  characteristic  of  finite  arrays. 
Unfortunately,  their  strength  is  also  the  source  of  their  weakness.  Since  the  current  on  each  array 
element  in  the  finite  dimension  is  an  unknown  quantity,  these  methods  are  limited  to  relatively  small 
arrays  by  available  computer  resources.  Many  authors  have  proposed  various  schemes  to  alleviate 
this  fundamental  problem  [4,  6, 7, 15, 22,  36, 37],  but  the  literature  is  devoid  of  calculated  results  for 
finite  arrays  larger  than  tens  or  hundreds  of  elements. 

The  limitations  of  these  finite  and  infinite  array  methods  suggest  a  hybrid  approach  for  the  large 
array  problem.  Surprisingly,  only  a  small  number  of  papers  on  such  techniques  have  been  published. 
Researchers  at  the  University  of  Illinois  [8]  have  investigated  truncated  and  curved  strip  arrays  via  a 
technique  which  approximates  the  induced  current  on  each  element  as  if  it  were  a  member  of  an  infinite 
array.  The  authors  then  replace  the  edge  element  currents  by  those  calculated  for  a  much  smaller 
array.  Unfortunately,  this  method  gives  poor  results  for  incidence  angles  near  grazing  where  the 
coupling  between  the  edge  elements  and  those  interior  to  the  array  becomes  significant.  In  addition, 
the  authors  make  no  attempt  to  model  the  dielectric  coatings  common  in  radome  applications. 

Felsen  and  Garin  have  recently  investigated  finite  wire  arrays  using  a  Floquet  current  approx¬ 
imation  on  two  semi-infinite^  sub-arrays  [10,  12].  They  combined  the  two  sub-arrays  such  that 

^The  term  “semi-infinite,”  in  the  context  of  this  dissertation,  will  always  refer  to  an  array  with  only  one  edge  analogous 
to  a  half-plane,  {i.e.  Along  the  semi-infinite  axis,  the  elements  go  on  “forever”  in  one  direction,  but  terminate  at  some 
point  in  the  opposite  direction.)  This  definition  contrasts  with  one  commonly  used  in  the  literature  to  refer  to  an  array 
which  has  a  finite  number  of  elements  in  one  direction  and  an  infinite  number  of  elements  in  the  orthogonal  direction. 
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the  composite  currents  are  that  of  the  desired  finite  wire  array.  Several  assumptions  constrain  their 
analysis.  Namely,  they  assume  all  propagation  is  ray  optical.  Also  they  ignore  any  deviation  in  the 
edge  element  currents  as  well  as  multiple  scattering  between  array  elements.  As  a  result,  their  for¬ 
mulation  gives  poor  results  near  the  grazing  incident  angles  common  for  forward  aspect  illumination 
of  radomes. 

Finally,  in  his  1973  paper,  Wasylkiwskyj  [40]  examines  the  mutual  coupling  effects  in  semi¬ 
infinite  antenna  arrays.  His  analysis  is  based  on  a  Wiener-Hopf  solution  technique  for  an  infinite  order 
difference  equation  modeling  the  antenna  element  port  currents.  Wasylkiwskyj ’s  transform  approach 
provides  an  alternate  solution  methodology  which,  to  the  author’s  knowledge,  has  not  been  applied 
to  large  finite  arrays. 

The  need  to  account  for  array  truncation  in  large  arrays  especially  at  grazing  incident  angles 
along  with  the  short-comings  of  existing  hybrid  finite-infinite  analysis  techniques  motivates  this 
present  research.  Specifically,  by  more  accurately  modeling  the  currents  induced  on  a  semi-infinite 
array,  the  scattering  from  an  arbitrarily  large  finite  array  can  be  synthesized  in  a  manner  analogous  to 
that  used  by  Felsen  and  Carin.  The  resulting  formulation  should  provide  more  accurate  radar  cross 
section  (RCS)  predictions,  especially  in  the  areas  where  existing  hybrid  techniques  fail. 

1 . 1  Problem  Statement 

Improved  FSS  analysis  methods  are  required  for  the  design  of  advanced  antenna  radomes 
used  in  modem  aircraft.  Existing  techniques  based  on  infinite  array  models  cannot  account  for 
the  perturbations  caused  by  edge  effects  at  radome-to-aircraft  boundaries  nor  the  radome  curvature 
characteristic  of  conformal  designs.  Finite  array  analysis  methods,  while  able  to  handle  these 
phenomena,  are  restricted  to  relatively  small  array  sizes  by  modem  day  computer  memory  limitations, 
especially  if  dielectric  layers  are  involved.  The  infinite/finite  hybrid  techniques  proposed  to  date  may 
not  sufficiently  represent  the  tme  coupling  between  the  edge  elements  and  the  central  portion  of  a 
large  array.  The  objectives  of  this  research  are  first,  to  develop  and  test  a  new  hybrid  PMM/MM 
approach  which  addresses  this  coupling  issue  to  more  accurately  predict  the  performance  of  large 
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finite  planar  arrays  and  second,  to  use  this  approach  to  analyze  the  isolated  “diffraction”  of  a  single 
FSS  termination. 

1 .2  Solution  Approach 

This  dissertation  presents  a  new  approach  to  combining  infinite  array  and  small  finite  array 
analysis  based  on  the  concept  of  a  “physical  basis  function”  (PBF).  Although  couched  in  different 
terms,  the  PBF  concept  was  originally  proposed  by  Morita  [24]  and  later  Tew  and  Tsai  [38]  in  an 
effort  to  apply  the  MM  to  infinite  planar  structures.  Burnside,  Yu,  and  Marhefka  [5]  then  applied  it 
to  the  wedge  diffraction  problem,  combining  the  geometrical  theory  of  diffraction  (GTD)  with  the 
MM.  To  the  author’s  knowledge,  however,  the  present  work  represents  the  first  application  of  the 
PBF  concept  to  periodic  structures. 

The  core  of  the  PBF  concept  involves  a  fairly  straight-forward  application  of  a  priori  knowledge 
of  the  induced  current  behavior  to  the  MM.  Consider  the  currents  impressed  on  a  large  planar  array 
by  a  plane  wave  illumination  (large  implying  many  wavelengths).  Starting  from  an  edge  and  moving 
toward  the  center  of  the  array,  there  will  be  some  point  where  the  current  amplitudes  become 
approximately  constant  over  a  large  central  section  of  the  finite  array.  Over  this  portion  one  can  use 
the  PMM  to  determine  the  representative  admittance  or  impedance  of  these  elements.  This  PBF  can 
then  be  used  along  with  the  edge  element  basis  functions  to  create  a  coupling  matrix  which  relates 
the  unknown  current  expansion  coefficients  to  known  coefficients  dependent  on  the  incident  field. 
The  coupling  matrix  now  explicitly  incorporates  the  influence  of  the  “infinite”  portion  of  a  large 
finite  array  into  the  small  finite  array  analysis  in  a  manner  analogous  to  how  the  interconnection 
matrix  of  diakoptic  theory  [13]  combines  the  impedances  of  individual  structural  elements  into  a 
composite  structure  matrix.  This  technique  represents  a  departure  from  traditional  hybrid  methods 
which  calculate  the  currents  mduced  on  the  two  portions  of  the  array  separately  and  then  combine  their 
effects  in  the  far-field  through  a  radiation  integral  analogous  to  the  physical  optics  approximation. 
The  remainder  of  this  section  lays  out  the  general  hybrid  PMM/MM  approach  using  a  wire  array  in 
free-space  as  an  example  geometry. 
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Consider  the  semi-infinite  array  of  thin  wires  shown  in  Figure  2a.  Note  that  the  wires  in  the 
array  are  infinitely  long  and  form  an  array  which  continues  forever  to  the  right.  Using  the  surface 
equivalence  theorem,  one  defines  a  new  geometry  where  impressed  currents  radiate  into  free-space. 
Specifically,  one  replaces  the  perfect  electric  conductor  (PEC)  with  a  Huygens’  surface  and  defines 
impressed  surface  currents  of  the  form 


Js  =  hxH 
M  ^  =  Exh 


(1) 


where  and  Mg  are  the  impressed  electric  and  magnetic  surface  currents,  n  is  the  outward  normal 
of  the  Huygens’  surface,  and  H  and  E  are  the  magnetic  and  electric  fields  at  the  Huygens’  surface 
in  the  original  problem.  With  a  thin  wire  approximation,  the  electric  surface  currents  are  modeled  by 
z  directed  filamentary  currents  at  the  center  of  the  wires.  There  are  no  impressed  magnetic  currents 
since  the  total  tangential  electric  field  is  identically  zero  on  the  PEC  surface  in  the  original  problem. 
The  resulting  equivalent  problem  is  depicted  in  Figure  2b. 

The  equivalent  problem  contains  an  infinite  number  of  unknown  currents  making  a  numerical 
solution  via  a  traditional  MM  approach  impossible.  In  order  to  make  the  problem  tractable,  one 
assumes  the  current  amplitudes  vary  significantly  only  over  the  first  B  elements.  The  remaining 
elements  have  approximately  equal  amplitude  currents  which  differ  by  a  Floquet  phasing  factor 
dependent  on  the  incident  plane  wave  direction.  This  assumption  reduces  the  problem  to  that  of 
finding  (B  -|- 1)  unknowns  as  shown  in  Figure  3. 

One  next  forms  an  integral  equation  relating  the  unknown  currents  to  the  known  incident  field. 
As  is  common  in  these  formulations,  the  total  tangential  electric  field  boundary  condition  is  used  to 
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Figure  2.  Semi-Infinite  Wire  Array 
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Figure  3.  Equivalent  Geometry 


obtain  an  integral  equation  of  the  form 


where  n  is  the  outward  normal  of  the  Huygens’  surface,  E  is  the  incident  electric  field  on  the 
wire  surfaces  {S),  and  E^{Jb)  is  the  scattered  electric  field  due  to  the  6*^  current  element.  Note 
that  in  general,  one  can  directly  express  the  fields  due  to  a  current  via  the  radiation  integral  and  the 
appropriate  dyadic  Green’s  function  [39] 
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where  r{R,E!)  is  the  dyadic  Green’s  function,  Jb{I^)  is  the  6*^  element  current  density,  and  r! 
and  R  denote  the  source  and  observation  locations  respectively. 

One  now  solves  Equation  (2)  via  the  MM  by  first  expanding  the  unknown  currents  with 
appropriate  known  basis  functions,  4>b{p)-  In  this  case,  B  impulse  basis  functions  and  one  PBF  are 
required.  These  take  the  form 

Jb{p')  «  JbMp')  (4) 


where 


and 


'  oo 

_  ~  b  =  0 

Mp  )  =  z  <  9=0 

^  S{p'-p,)  6=  1,2,. ..5 

Pq  =  x  qd^ 


(5) 


(6) 


Note  that  locates  the  “edge”  of  the  PBF  and  s  denotes  the  incident  plane  wave  propagation 
direction.  Inserting  Equation  (4)  in  Equation  (2),  one  obtains  the  expression 


B 

ny.W  =  E\Jb'^i,{p')) +  e[p,J^,Ji,...JB)  (7) 

b=0  p^s 

where  e  is  the  error  function  associated  with  the  current  approximation. 

Following  Harrington’s  MM  formulation  [18],  the  next  step  in  the  numerical  solution  involves 
defining  an  inner  product  and  forcing  a  chosen  set  of  (5 + 1 )  testing  functions,  Ot  (p) ,  to  be  orthogonal 
to  the  error  function.  For  the  wire  geometry,  the  filamentary  current  nature  suggests  a  Dirac  delta 
testing  function.  In  order  to  avoid  the  Green’s  function  singularity,  one  places  the  test  location  a 
wire  radius  away  from  the  basis  function  for  the  self-coupling  terms.  One  also  defines  the  testing 
function  associated  with  the  PBF  to  be  nonzero  only  over  a  chosen  reference  wire  element  to  keep 
the  coupling  matrix  elements  bounded.  Dotting  these  testing  functions  into  Equation  (7)  and  taking 
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the  inner  product  results  in  the  expression 


J  nxF(p)]  •0t(p)(ip  =  -J 


B 


hxY,E‘(JiMp)) 


s 

B 


6=0 


et{p)  dp 


=  -E^‘/ 

6=0  i 


(8) 


hxE‘’{(j>^{p))  \  ■9t{p)dp 


where  t  =  0, 1, . . .  B  and  E"{(f>^^[p))  is  the  field  at  p  radiated  by  the  basis  function.  Equation  (8) 
now  defines  a  (5  +  1)  x  (B  +  1)  matrix  equation  which  can  be  expressed  in  matrix  form  as 


Zoo 

Zob 

Jo 

.  . 

Zto 

Ztb 

.  . 

where  the  0  subscript  denotes  the  basis  and  testing  modes  associated  with  the  PBF.  Clearly  the 
“influence”  of  the  neighboring  infinite  array  on  the  edge  elements  is  contained  in  the  cross-coupling 
terms,  Zo6  and  ZtQ.  Note  the  terms  with  the  t  and  b  subscripts  actually  represent  sub-arrays  and  sub¬ 
matrices  associated  with  the  edge  elements.  One  can  calculate  the  entries  in  the  matrix  Equation  (9) 
for  the  wire  array  from  Equation  (3)  by  using  the  appropriate  dyadic  Green’s  function. 


Vo  = 

e-j^^Po 

(10) 

Vt  = 

(11) 

o 

II 

OO 

(12) 

g=0 

II 

(13) 

Zto  — 

OO 

(14) 

q=0 

Ztb  = 

(15) 
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One  can  also  invoke  the  reciprocity  theorem  [17]  to  exploit  symmetries  in  the  edge  element  sub-matrix 
thereby  reducing  the  matrix  fill  time.  Standard  matrix  equation  techniques  can  then  be  employed  to 
solve  the  system,  obtaining  the  desired  current  coefficients. 

The  final  step  in  the  hybrid  PMM/MM  procedure  involves  using  the  calculated  current  ex¬ 
pansion  coefficients  to  determine  the  scattered  fields  via  a  radiation  integral.  Specifically,  one  finds 
the  field  from  a  large  finite  array  by  solving  the  semi-infinite  array  problem  twice,  once  for  each 
edge.  One  then  synthesizes  the  currents  on  the  finite  array  by  shifting  the  edge  element  currents 
appropriately  and  assigning  the  “central”  elements  the  PBF  current  amplitude  with  an  appropriate 
linear  phase.  The  radar  echo  width  is  then  calculated  from  the  scattered  fields. 

1.3  Dissertation  Summary 

The  remaining  chapters  take  the  general  procedure  just  outlined  and  apply  it  to  two  geometries; 
a  two-dimensional  wire  array  and  a  three-dimensional  slot  array,  both  in  a  stratified  dielectric  media. 
Chapter  II  begins  by  deriving  the  coupling  matrix  entries  for  the  wire  array.  It  presents  both  spatial 
and  spectral  domain  forms  for  the  radiation  from  the  PBF  and  edge  element  basis  functions  as  well  as 
addresses  the  complexities  introduced  by  adding  a  dielectric  slab  to  the  geometry.  Next,  the  chapter 
formulates  the  excitation  vector  entries,  again  taking  into  account  the  slab  boundaries.  Finally,  the 
chapter  concludes  with  a  discussion  of  how  the  solution  vector  entries  can  be  used  to  calculate  a 
far-field  radiation  pattern  for  a  large  finite  array.  Chapter  III  gives  the  corresponding  analysis  for 
the  slot  array.  Chapter  IV  examines  the  validity  of  the  PBF  formulation  through  reference  code 
and  measurement  comparisons.  In  general,  the  array’s  induced  currents  are  compared,  since  radar 
echo  width  (REW)  calculations  tend  to  smooth  out  errors  in  the  solution  vector.  For  the  measurement 
comparisons,  however,  the  author  was  forced  to  use  RCS  calculations  as  only  measured  RCS  data  was 
available.  Finally,  Chapter  V  re-caps  the  dissertation,  summarizing  the  salient  results  and  suggesting 
areas  for  future  research. 
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IL  Two-Dimensional  Semi-Infinite  Wire  Array 


A  journey  of  a  thousand  miles  begins  with  a  single  step. 

-  A  Chinese  Proverb 

As  a  proof-of-concept,  this  chapter  develops  the  hybrid  PMM/MM  solution  given  by  Equa¬ 
tion  (9)  for  the  simple  two-dimensional  wire  geometry  depicted  in  Figure  4.  The  resulting  expressions 
form  the  basis  for  the  wire  computer  code  validated  in  Chapter  IV.  Equation  (9)  naturally  breaks 
down  into  three  parts.  This  chapter  follows  this  organization  by  first  examining  the  coupling  matrix, 
followed  by  the  excitation  vector,  and  finally  the  solution  vector.  Where  appropriate,  each  of  these 
sections  is  divided  into  subsections  covering  the  homogeneous  and  stratified  media  cases  respectively. 

2.1  Coupling  Matrix 

The  free-space  analysis  used  to  illustrate  the  PBF  concept  in  Chapter  I  forms  the  common 
starting  point  for  both  the  dielectric  coated  wire  and  slot  arrays.  In  the  wire  case,  the  array  is  in  a 


Figure  4.  Semi-Infinite  Wire  Array  Geometry 
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homogeneous  media  with  region  two’s  constitutive  parameters  (see  Figure  4).  Since  the  filamentary 
currents  are  z  directed,  one  can  write  Equation  (8)  as  the  scalar  equation 


T 


oo  ^ _ 

.  «=o 


6=1 


(16) 


where  'p^  (test  function  location)  is  confined  to  the  Huygens  surface  of  the  wire,  (PBF  element 
location)  is  confined  to  the  center  of  the  5*^  wire,  pj,  (basis  function  location)  is  confined  to  the  center 
of  the  wire,  and  all  quantities  related  to  the  constitutive  parameters  are  those  of  region  two.  Note 
how  the  shifting  property  of  the  Dirac  delta  function  reduces  the  surface  integrals  to  point  evaluations. 
Equation  (16)  forces  equality  between  the  fields  radiated  by  the  equivalent  currents  and  the  incident 
field  at  these  points  on  the  Huygens  surface. 

Equation  (16)  is  known  as  a  “spatial  domain”  expression.  In  order  to  account  for  the  dielectric 
slab,  one  could  derive  a  non-ffee-space  Green’s  function  satisfying  the  appropriate  field  boundary 
conditions  at  the  slab  boundaries.  The  author  chose  an  alternative  approach  based  on  a  “spectral 
domain”  representation  of  Equation  (16).  In  the  spectral  domain,  one  deals  with  plane  wave  field 
expansions.  Thus,  planar  boundaries  can  be  introduced  with  relative  ease  via  Fresnel  coefficients  [17]. 
In  addition,  the  spectral  domain  summations  are,  in  general,  more  convergent  than  their  spatial  domain 
counterparts,  leading  to  more  efficient  computer  codes.  The  means  for  transforming  to  the  spectral 
domain  vary  for  the  different  parts  of  Equation  (16).  Specifically,  one  can  partition  the  right  hand 
side  into  two  parts  representing  the  fields  radiated  by  the  PBF  and  the  fields  radiated  by  the  individual 
edge  elements.  The  following  subsections  develop  the  spectral  domain  representation  of  these  fields 
which  are  related  to  the  coupling  matrix  entries  through  the  unknown  current  coefficients,  Jj. 


2.1.1  Physical  Basis  Function. 

2. 1.1.1  Physical  Basis  F unction  Radiating  to  Test  Function  in  a  Homogeneous  Media. 
Figure  5  depicts  the  PBF  in  a  homogeneous  media  with  constitutive  parameters  and  /i2.  One  can 
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Figure  5.  PBF  Radiating  to  Test  Function  Location  in  Homogeneous  Media 


express  the  fields  radiated  by  the  PBF  in  the  spatial  domain  as  a  collection  of  cylindrical  waves. 

j  OO  _ 

g=0 

=  E  + !/?)  e-’'’*”''-  (17) 

Taking  a  hint  from  Munk’s  work  [26,  28],  one  would  like  to  use  the  Poisson  sum  formula  to  convert 
Equation  (17)  to  the  spectral  domain.  Unfortunately,  the  one-sided  summation  prevents  a  direct 
application  of  the  Poisson  sum  formula.  Appendix  A  presents  a  one-sided  version  of  the  regular 
Poisson  sum  formula  [29]  applicable  to  summations  such  as  in  Equation  (17).  Equipped  with  this 
formula,  one  can  convert  Equation  (17)  to  a  collection  of  plane  waves.  Specifically,  one  uses  the 
transform  pair  [3] 
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Note  that  both  Py  and  Ty  are  defined  to  be  either  positive  real  or  negative  imaginary  in  order  to  satisfy 
the  radiation  boundary  condition,  (i.e.  The  fields  must  be  bounded  at  infinity.)  This  also  restricts  yt 
to  the  upper  half-plane. 

Equation  (21)  gives  the  spectral  domain  representation  of  the  fields  radiated  by  the  PBF  in  a 
homogeneous  media.  As  promised,  this  form  is  more  convergent  than  the  spatial  domain  form  owing 
to  the  exponential  damping  in  the  higher  order  spectral  terms.  Specifically,  when  -f-  >  1  and 

Ira; I  >  1  the  plane  waves  become  evanescent,  decaying  rapidly  for  test  locations  out  of  the  array 
plane.  It  is  also  interesting  to  note  the  one-sided  Poisson  sum  generates  both  discrete  and  continuous 
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spectra.  The  discrete  spectrum  incorporates  the  periodic  nature  of  the  array  whereas  the  continuous 
spectrum  accounts  for  the  diffraction  effects  caused  by  the  array  termination.  In  fact,  the  similarity 
of  the  continuous  spectral  integral  to  the  diffraction  integral  encountered  in  the  half-plane  scattering 
problem  provides  clues  to  its  efficient  evaluation. 

The  first  of  these  clues  allows  one  to  remove  the  weak  singularity,  from  the  integrand. 
Felsen  and  Marcuvitz  [11]  suggest  the  variable  substitution 


Tx 

=  sin(z/) 

(24) 

=  cos(i/) 

(25) 

dTx 

=  cos(i/)  dTx 

(26) 

leading  to  the  angular  spectral  integral 


Ik 


-i 


sin(i/)+V(  cos(i/)] 

(sin(i/)  -  p^) 


dv 


(27) 


The  integration  contour,  C,  lies  in  the  regions  of  the  complex  u  plane  where  Im{cos{u)}  <  0  as 
indicated  by  the  shaded  areas  in  Figure  6.  The  deviation  of  C  from  the  real  axis  represents  a  slight 
loss  in  the  wave  number  /?.  Note  the  pole  at  Ug  =  sm~^{px)  is  not  shown  in  the  figure.  One  can 
further  simplify  Equation  (27)  by  defining  a  new  cylindrical  coordinate  system  displaced  from  the 
origin  by  one  half  the  inter-element  spacing  (Figure  7). 

With  the  definitions 


Pt 


tan  ^ 


(28) 

(29) 
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Im(  V) 


Figure  6.  Angular  Spectral  Integral  Contour 


and  the  cosine  addition  formula,  the  angular  spectral  integral  becomes 


h 


/ 


Q-j^Pt  coa{v-it) 
(sin(i/)  -  p^) 


du 


(30) 


One  is  now  faced  with  the  task  of  evaluating  Equation  (30).  Even  though  it  is  a  principle  value 
integral,  the  non-integrable  singularity  at  Vg  makes  straight-forward  numerical  integration  difficult. 
Again  taking  a  clue  from  its  similarity  to  a  diffraction  integral,  one  might  consider  evaluating 
Equation  (30)  asymptotically.  In  fact,  for  many  of  the  off-diagonal  coupling  matrix  terms,  /5pt 
can  be  considered  a  large  parameter.  Appendix  C  presents  such  an  asymptotic  evaluation,  clearly 
demonstrating  the  tremendous  computational  advantage  over  numerical  evaluation  schemes. 

Where  asymptotic  evaluation  is  not  possible,  one  must  evaluate  Equation  (30)  numerically. 
The  author  found  the  following  singularity  extraction  technique  works  quite  well.  One  first  breaks 
the  integration  contour  into  three  regions,  the  vertical  path  [— (|  -|-  joo),  —  |),  the  real  axis  (-f ,  f ), 
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Figure  7.  Alternate  Cylindrical  Coordinate  System 


and  the  vertical  path  (|,  |  +  joo].  Along  the  vertical  paths,  one  redefines  Uoid  =  T  (f  +  J^new) 
where  0  <  Vnew  <  oo,  transforming  them  to  the  real  axis.  Thus,  the  angular  spectrum  transformation 
for  these  paths  becomes 

=  T  sin  cosh(z/)  (31) 

Equation  (30)  now  can  be  written 


=  -^1 


giPpt  siii(ji/+0<) 

(cosh(i^)  +  px) 


g-j0Pt  cos{v-<ft) 

(sin(i/)  -  p^) 


du  +  j- 


(cosh(z/)  -  px) 


h  +  h  +  h 


where  the  principle  value  occurs  in  only  one  integral  for  a  particular  value  of  px. 

One  next  forces  the  numerator  of  the  appropriate  singular  integral  to  go  to  zero  exactly  at  the 
singular  point,  effectively  removing  the  singularity.  The  extracted  singular  integral  is  then  in  a  form 
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suitable  for  an  analytic  solution.  Thus,  Equation  (32)  becomes 


h  = 


- - - — r - r -  du  -|-  N\q~T  7 - T  T 

(cosh(^')  +  px)  J  (cosh(z/)  +  px) 

Lo  0 


00 

/ 


+ 


+ 


f(‘ 


I 

“■  2 
00 

/ 


-jPpt  cos(i/-(^i) 


(sin(i/)  -  px) 

(cosh(z^)  -  Px) 


-  N2o) 

^du  +  N2 


4 


dv 


2 

OO 


dv  +  N?, 


(sin(z/)  -  px) 
dv 

(cosh(z/)  -  px) 


(33) 


where 


iVlo 

=  —j  ^j^Pt  I'o + '?^< ) 

(34) 

iV2„ 

—  g-jfiptCoa{vo-4t) 

(35) 

iV3„ 

—  j^i^Pi  ein{jvo-^t) 

(36) 

Abramowitz  and  Stegun  [1]  provide  analytic  solutions  for  the  extracted  integrals 


OO 

/ 


dv 


(cosh(z/)  +  Px) 


2 

/ 


dv 


(sin(i/)  -  Px) 


OO 

/ 


dv 

(cosh(i/)  -  Px) 


tan 


-1 


Pi 


1-f />X 


0 

2 

2 


In 


l+^r  +  \/ 
I+Pi-^/pE-I 


-  tan-^  f 


, _ tan-M^^ 

l-/>r  +  -v/p|^ 


In 


(37) 


(38) 


(39) 
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where  one  chooses  either  the  upper  or  lower  line  for  <  1  or  p^  >  1,  respectively.  Thus,  the 
spectral  domain  expression  for  the  fields  radiated  in  a  homogeneous  media  can  be  written 


rpPBFf-^  N  _  VJo 


E 

k=^oo 


j  qJPPx  -f" 

+  - - h 

TT 


(40) 


where  h  is  evaluated  either  asymptotically  or  numerically. 

Equations  (17)  and  (40)  provide  alternative  representations  for  the  field  radiated  by  the  PBE 
These  expressions,  with  the  Jq  constant  removed,  define  the  mutual  coupling  between  the  basis  and 
testing  functions.  Since  filling  the  coupling  matrix  constitutes  a  significant  portion  of  the  overall 
computation  time,  one  would  like  to  use  the  most  efficient  form.  The  next  few  paragraphs  examine 
both  the  spatial  and  spectral  domain  expressions  for  the  free-space  semi-infinite  wire  array  Green’s 
function  calculated  from  Equations  (17)  and  (40)  by  multiplying  with  the  constant  Figure  5’s 
geometry,  with  parameters  =  O.lAo,  s  =  —[x  cos(30“)  -|-  y  sin(30°)],  and  (€2,^2)  =  (^0,^0), 
serves  as  a  typical  test  case.  The  test  location,  p^,  is  scanned  along  a  path  parallel  to  the  x  axis  at  an 
“elevation”  of  yt  =  O.lAo. 

Figure  8  contains  the  Green’s  function’s  real  and  imaginary  components  and  magnitude  as  the 
test  point  sweeps  past  the  array  edge.  As  one  would  expect,  the  magnitude  of  the  Green’s  function 
decays  monotonically  to  the  left,  away  from  the  array  edge.  Conversely,  as  the  test  point  moves  to 
the  right  over  the  array,  the  magnitude  settles  down  to  the  same  value  one  obtains  from  a  fully  infinite 
array.  Note  that  the  data  in  this  figure  can  come  from  either  the  spatial  or  spectral  domain  expressions. 
The  question  is  which  one  is  more  efficient. 

Figure  9  shows  the  number  of  terms  required  by  each  expression  to  converge  to  the  criteria 


Sum  —  Truth 
Truth 


<  0.0001 


(41) 


The  “Truth”  represents  the  value  obtained  from  either  expression  pushed  to  the  point  where  the  next 
term  does  not  alter  the  previous  partial  sum  by  more  than  machine  accuracy.  In  addition,  the  author 
used  the  “circle-average”  summation  acceleration  technique  described  in  Appendix  B  to  reduce  the 
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Figure  8.  Semi-Infinite  Wire  Array  Green’s  Function 

number  of  terms  required  by  the  spatial  domain  expression.  The  spectral  domain  summation  did 
not  require  acceleration.  Clearly,  the  spectral  domain  is  the  domain  of  choice  if  one  is  interested  in 
reducing  the  number  of  summation  terms. 

The  final  figure  in  this  section  numerically  confirms  the  proper  evaluation  of  the  continuous 
spectral  integral  in  Equation  (40).  From  Equations  (174)  and  (185)  in  Appendix  A,  it  is  clear  that,  as 
the  test  point  moves  to  the  right  over  the  semi-infinite  array,  the  continuous  spectrum  should  settle 
down  to  a  value  equal  to  the  discrete  spectrum  which,  in  turn,  is  exactly  half  the  value  obtained  for 
a  fully  infinite  array.  Figure  10  plots  the  discrete  and  continuous  spectra  in  the  complex  plane  as  the 
test  point  moves  on  the  scan  path  from  l.OAo  to  the  left  of  the  array  edge  to  I.OAq  to  the  right  of  the 
array  edge.  The  points  designated  by  the  “+”  show  the  discrete  and  continuous  spectra  have  nearly 
the  same  magnitude  and  phase  over  the  array.  As  the  test  point  moves  past  the  array  edge  (“x”)  to  a 
location  left  of  the  array  (“o”),  notice  how  the  magnitudes  again  become  similar,  but  now  are  strongly 
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X  position  (lambda) 

Figure  9.  Number  of  Terms  Required  for  Convergence 

“out  of  phase.”  The  resulting  cancelation  reduces  the  field  magnitude  as  the  test  point  moves  away 
from  the  array. 

To  re-cap,  this  section  first  developed  and  then  examined  the  spatial  and  spectral  domain 
expressions  for  the  mutual  coupling  terms  associated  with  the  PDF  in  a  homogeneous  media.  The 
results  clearly  demonstrate  that,  for  most  cases  where  the  test  location  is  out  of  the  plane  of  the  array, 
the  spectral  domain  expression  is  preferred.  The  next  section  addresses  the  geometry  containing  a 
dielectric  slab.  There  the  other  advantage  of  the  spectral  domain  form  will  be  apparent,  namely,  the 
ease  with  which  one  can  introduce  planar  boundaries. 

2.1.12  Physical  Basis  Function  Radiating  to  Test  Function  in  a  Stratified  Media.  As 
mentioned  earlier,  the  homogeneous  media  analysis  forms  the  foundation  for  introducing  a  layered 
structure.  Consider  adding  a  dielectric  slab  to  Figure  5.  In  addition  to  the  direct  coupling  between  the 
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Figure  10.  Complex  Spectra 

PBF  and  the  test  function  at  p^,  one  must  now  account  for  the  plane  waves  reflected  back  and  forth 
between  the  two  boundaries  (Figure  11).  Fortunately,  these  boundaries  are  infinite  parallel  planes. 
One  can  use  the  concept  of  effective  Fresnel  coefficients  [17]  to  calculate  the  plane  waves’  magnitude 
and  phase  as  they  reflect  off  the  boundaries.  In  essence,  the  problem  comes  down  to  determining  the 
composite  effect  of  a  plane  wave  radiated  by  the  PBF  rattling  between  the  slab  boundaries  and  then 
terminating  at  the  test  location. 

In  his  work  on  periodic  surfaces  in  stratified  dielectric  media,  Munk  [26,  28]  developed  a  set 
of  effective  reflection  and  transmission  coefficients  useful  for  propagating  a  plane  wave  through  an 
arbitrary  number  of  parallel  dielectric  slabs.  His  coefficients  use  a  polarization  decomposition  based 
on  the  plane  of  incidence  shown  in  Figure  12.  According  to  Munk’s  convention,  r]  lies  in  the  plane  of 
incidence  and  is  measured  from  the  negative  y  axis,  whereas  a  lies  in  the  xz  plane  and  is  measured 
from  the  x  axis.  The  author  adopts  Munk’s  convention  for  this  research.  The  interested  reader  is 
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Figure  11.  PBF  Radiating  to  Test  Function  Location  in  Stratified  Media 


referred  to  the  cited  reports  for  details  concerning  the  coefficients’  development.  Since  the  incident 
plane  wave  propagation  vector  is  confined  to  the  xy  plane,  the  z  polarized  fields  are  perpendicularly 
polarized  and  require  the  following  reflection  and  transmission  coefficients^ 


'^1,1+1 


V{l+l)Qy{l)  +  VinQyil+l) 


(42) 


r,. 


/+! 


^  'n(l+l)Qy(l)  -  V(l)gy{l+1) 
riH+i)Qy(i)  +  r]{i)Qy{i+i) 


(43) 


where  is  the  plane  wave  propagation  vector’s  y  component  in  the  slab  and  r]^i)  is  the 
slab’s  impedance.  Note  these  coefficients  apply  to  half-space  boundaries  only.  For  the  multi-layer 
geometries  discussed  in  Chapter  III,  one  must  use  the  effective  Fresnel  coefficients  defined  by  Munk. 


^  These  are  valid  for  waves  propagating  in  a  direction  with  a  positive  y  component.  One  can  calculate  the  reflection 
coefficient  for  waves  with  a  negative  y  component  from  the  relationship,  Fj+j,;  =  -Fi^j+i. 
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Figure  12.  Plane  Wave  Polarization  Decomposition 


Now,  to  determine  the  composite  plane  wave  at  one  groups  the  bouncing  plane  waves  into 
the  four  infinite  summations  shown  in  Figure  13.  These  summations  are  defined  in  terms  of  the  plane 
wave’s  departure  and  arrival  directions.  For  example,  summation  one  consists  of  the  plane  waves 
leaving  the  PBF  reference  location,  (xo,yo).  with  a  positive  y  component  and  arriving  at  the  test 
location,  {xt,  yt),  with  a  positive  y  component.  Summation  two  consists  of  the  plane  waves  with  a 
negative  y  departure  and  positive  y  arrival.  Summation  three  consists  of  plane  waves  with  a  positive 
y  departure  and  a  negative  y  arrival.  Finally,  summation  four  consists  of  the  plane  waves  with  a 
negative  y  departure  and  a  negative  y  arrival.  Mathematically,  these  summations  can  be  expressed 


^direct  +  .  .  .] 

^direct  \j^^^Q-iP2'2{.yh-yi)Qy2  _j_  +  ] 

^direct  ]j^^^Q-jl322{y2-yt)Qy2  _j_  p2^p2jg-i/?22(2v2-i/i-?/t)ej2  _(_  _  _ 

^direct  j■p2J^P2gg-J^22(2/6-^/^+2/2-3/l)e!(2  p2^p2^g-J^22(3/i,-i/t+2(2/2-l/l))e!(2  _|_ 


(44) 

(45) 

(46) 

(47) 
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1st  Summation  2nd  Summation  3rd  Summation  4th  Summation 


Figure  13.  Decomposition  of  the  Plane  Waves  Internally  Reflected  in  the  Dielectric  Slab 

where  is  the  plane  wave  propagating  directly  to  the  test  location  without  any  bounces.  In 

particular,  is  a  spectral  component  from  either  the  discrete  or  continuous  spectra  of  the 

homogeneous  media  Equation  (21).  The  author  found  it  numerically  expedient  to  subtract  out  the 
direct  term  from  Equation  (44)  and  express  it  in  the  spatial  domain.  The  remaining  summations  are  in 
the  form  of  geometric  series  which  can  be  expressed  in  closed  form.  Combining  these  summations, 
one  arrives  at  a  multiplier  to  the  homogeneous  media  field  expression  which  accounts  for  the  presence 
of  the  stratified  media.  Munk  calls  this  multiplier  a  “T-factor.”[26,  28] 
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One  can  now  express  the  field  radiated  by  the  PBF  in  the  dielectric  slab  as 


Er^iPt) 


/?2^2>^0 


Lq=0 


g  3p2p2'Pt 


E{py2) 


+ 


P2dx  I  Py2 

j  '7  Q-iP2'f'-P^Q-3P2^{Tx-px2) 


'^yip'x  Px2^ 


T{ry)drx 


(48) 


where  the  T-factor^  is 


T{Qy2) 


r23e"-'^^2(3/2-yt)ea2  _  ri2e^^^23/iej,2)  _  +  P23e“^’^^22/2ey2^ 

1  +  ri2r23e-^2^2(3'2-3/i)e!,2 


(49) 


Since  Equation  (48)  is  based  on  the  homogeneous  expression,  the  spectral  integral  has  the  same  singu¬ 
larities  as  before.  One  can  use  the  same  regularization  procedures  to  remove  the  weak  ^  singularity 
and  the  singularities  associated  with  the  k  dependent  grating  lobes.  Unlike  the  homogeneous  media 
case  however,  one  must  deal  with  the  singularities  introduced  by  the  T-factor.  These  singularities 
occur  when  the  dielectric  slab  has  a  greater  permittivity®  than  the  surrounding  material,  causing  the 
spectral  terms  to  excite  surface  waves  in  the  dielectric  slab.  Quantitatively,  this  occurs  when  the 
denominator  in  Equation  (49)  goes  to  zero.  Lets  examine  this  situation  in  more  detail. 

Using  Equation  (23)  in  Equation  (49),  one  can  write  the  condition  for  a  surface  wave  singularity 
as 

1  +  =  0  (50) 

^The  expression,  as  written,  used  the  fact  the  PBF  is  located  at  yo  =  0. 

^Assuming  non-magnetic  materials. 
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where  from  Equation  (43) 


and  Tx2  sweeps  through  the  range  —  oo  <  Tx2  <  oo.  Recalling  that  Equation  (48)  is  based  on  a 
closed  form  representation  of  the  infinite  series,  one  might  be  concerned  with  its  validity  over  the 
entire  rx2  range.  Fortunately,  as  Skinner  [33]  indicated  in  his  dissertation,  one  can  add  a  slight  loss 
to  the  y  propagation  direction  and  thus  extend  the  equation’s  valid  range  to  the  entire  integration 
range.  This  approach  also  jives  with  the  physical  reality  that  no  material  is  perfectly  lossless.  For 
|r3.2|  <  min  y^)’  corresponding  plane  wave  forms  such  an  acute  angle  with  the  slab 

normal  that  some  of  the  wave  escapes  the  surface  (i.e.  |r|  <  1)  and  no  singularity  exists.  On  the  other 
hand,  when  |rj.2|  >  1,  the  y  propagation  becomes  purely  evanescent  and  the  reflection  coefficients 
can  be  written 


where  £2  >  €i,  63.  Thus, 


and  again  no  singularity  exists.  Therefore,  the  surface  wave  singularities  are  confined  to  the  range, 
min  ^y/^)  y^)  <  k®2|  <  1  and  can  be  found  by  solving  Equation  (50)  for  rx2- 


Since  Equation  (50)  is  a  transcendental  equation,  one  must  use  an  iterative  technique  such  as 
Newton-Raphson  to  locate  the  singularities.  Figures  14,  15,  and  16  contain  plots  of  the  T-factor  in 
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rhox 

Figure  14.  T-factor  Singularity  in  a  €2  =  2eo,  0.2Ao  Thick  Slab  with  a  0.1  Loss 

the  neighborhood  of  one  of  these  singularities.  The  geometry  in  these  plots  consists  of  a  0.2  Ao  thick 
slab  in  free-space  with  a  relative  permittivity  62  =  2eo  and  a  y  directed  loss  of  0.1,  0.01,  and  0.001 
respectively.  As  Skinner  [33]  found  in  his  work,  the  area  under  the  real  and  imaginary  curves 
approaches  a  limiting  value  as  the  loss  goes  to  zero.  Based  on  this  observation,  the  author  chose  to 
handle  the  T-factor  singularities  by  adding  a  slight  y  directed  loss  to  the  denominator  and  integrating 
the  continuous  spectrum  finely  in  the  neighborhood  of  the  former  singularity. 

Equation  (56)  gives  the  final  form  for  the  field  radiated  by  the  PBF.  The  author  implemented 
this  expression  in  the  two-dimensional  wire  code  validated  in  Chapter  IV. 
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Figure  15.  T-factor  Singularity  in  a  £3  =  2eo.  0.2Ao  Thick  Slab  with  a  0.01  Loss 


Nlo  =  (58) 

N2o  =  (59) 

iV3o  =  (60) 

As  before,  the  extracted  integrals  are  computed  using  Equations  (37),  (38),  and  (39)^.  Figure  17 
plots  the  PBF  magnitude  using  Equation  (56)  for  four  different  slab  permittivities  over  an  x  directed 
scan  similar  to  that  of  Figure  8.  In  this  case,  dx  =  0.3Ao,  S3  =  —[x  cos(30°)  +  y  sin{30°)],  and 

■^A  note  of  caution  concerning  the  singularity  extraction  procedure  for  the  grating  lobe  singularities;  The  numerator 
must  be  a  fairly  well-behaved  function  in  the  vicinity  of  the  singularity.  Thus,  in  order  to  use  the  extraction  procedure, 
one  must  ensure  the  T-factor  singularity  is  sufficiently  far  away  from  the  grating  lobe  singularity.  The  author  found  that, 
in  practice,  this  condition  is  not  unduly  restrictive. 
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Figure  17.  PBF  Radiation  in  the  Presence  of  a  Dielectric  Slab  of  Various  Relative  Permittivities 

yt  =  O.OSAo.  The  dielectric  slab  boundaries  are  at  yi  =  — O.lAo,  and  t/2  =  O.lAo  with  the  external 
regions  one  and  three  defined  to  be  free-space.  Finally,  a  0.01  T-factor  loss  tangent  ensures  the  field 
values  are  converged.  As  one  might  expect,  the  =  1  case  decays  monotonically  to  the  left  of  the 
array  edge.  Surface  waves  become  apparent  as  one  increases  the  relative  permittivity.  Also,  the  field 
curves  change  smoothly  as  the  relative  permittivity  increases,  giving  confidence  in  the  numerical 
implementation. 

The  final  figure  in  this  section  shows  the  effect  of  changing  the  “elevation”  of  the  scan  path 
above  the  array.  Figure  18  plots  the  PBF  fields  for  yt  =  O.OlAo,  0.025Ao,  and  O.OSAo  above  the  array 
in  a  free-space  dielectric  slab.  One  may  note  the  effect  of  the  individual  wires  on  the  field  as  the 
observation  point  approaches  the  array  plane.  As  expected,  the  field  magnitude  increases  dramatically 
in  the  neighborhood  of  each  wire.  Even  so,  the  field  values  remain  comfortably  bounded  with  yt 
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Figure  18.  PBF  Radiation  for  Various  Scan  Path  Elevations  above  the  Array  Plane 

as  close  as  O.OlSAo  away,  allowing  one  to  use  a  fairly  small  wire  radius  for  the  coupling  matrix, 
self-impedance  calculations. 

This  completes  the  PBF  subsection.  Equation  (56),  with  Jo  divided  out,  gives  the  impedance 
between  the  PBF  and  the  test  functions.  The  only  restriction  involves  keeping  the  test  location  out  of 
the  array  plane  for  the  self-impedance  terms.  The  subsection  addressed  issues  concerning  the  spectral 
integral’s  evaluation  along  with  the  difficulties  encountered  when  one  introduces  a  stratified  media. 
The  next  subsection  completes  the  coupling  matrix  section  with  the  development  of  an  expression  for 
the  radiation  from  an  individual  edge  element  in  a  stratified  media. 

2.1.2  Edge  Element  Basis  Function. 

2. 1.2.1  Edge  Element  Basis  Function  Radiating  to  Test  Function  in  a  Homogeneous 
Media.  Returning  to  Equation  (16),  one  can  express  the  field  radiated  by  a  single  edge  element  in 
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a  homogeneous  media  as 

-  pj)  (61) 

One  would  like  to  convert  Equation  (61)  to  the  spectral  domain  in  anticipation  of  introducing  a 
stratified  media.  In  order  to  accomplish  this,  the  author  took  what  may  first  appear  to  be  a  circuitous 
route.  Using  the  array  scanning  method  (ASM)  reported  by  Munk  and  Burrell[27],  the  author 
envisioned  a  “pseudo-array”  around  the  edge  element.  The  array  allows  one  to  introduce  the  Poisson 
sum  formula  to  transform  to  the  spectral  domain.  The  field  due  to  the  original  edge  element  can  then 
be  extracted  using  Fourier  methods  by  identifying  the  element  as  the  central  term  in  a  Fourier  series. 
Appendix  D  details  this  procedure. 

2. 1.2.2  Edge  Element  Basis  Function  Radiating  to  Test  Function  in  a  Stratified  Media. 

Based  on  the  homogeneous  media  development,  one  can  introduce  a  dielectric  slab  to  the  geometry 
in  exactly  the  same  manner  as  for  the  PBF.  In  fact,  the  T-factor  expression  developed  in  Section  2. i.i. 2 
can  be  used  without  modification.  As  before,  one  finds  it  expedient  to  decompose  the  edge  element 
field  into  direct  and  bounce  path  components 


56 ^  (p J  (62) 

with  kept  in  the  spatial  domain  and  expressed  as  a  spectral  summation. 

Equation  (61)  provides  the  direct  component  where  again  all  constitutive  parameters  are  those  of 
region  two.  One  can  use  the  ASM  to  derive  an  expression  for  the  bounce  component. 

One  first  defines  a  “pseudo-array”  around  the  edge  wire  located  at  pj, 

Er‘'(p,)  =  E  (63) 
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where 


Pqb  =  ®  (94  -{b-  1)4) 

(64) 

(65) 

In  order  to  isolate  the  surface  wave  pole  excitation  to  the  central  term  of  the  spectral  sum,  one  sets 
the  inter-element  spacing  to  Now,  from  Equation  (174)  in  Appendix  A  and  the  frequency 

shift  property  of  the  Fourier  transform,  one  can  transform  Equation  (63)  to  the  spectral  domain  form 


Ef'‘y(p,)  =  f; 


fc=— oo 


^ASM 


(66) 


where 


^ASM  ^  ^ 


„ASM  I 

^  /jt  I  ^ 

.  4 


(67) 


Using  the  spectral  form  in  the  ASM  integrand,  one  can  introduce  the  T-factor  to  account  for  the 
dielectric  slab.  Upon  interchanging  the  order  of  summation  and  integration,  the  final  expression  for 
the  bounce  term  in  Equation  (62)  can  be  written 


Tpboiince  (—  \  724 

2A7 


~  r 

S/ 


Hh  \p^^^ ■p^-p^^^{b~\)dx 


fyASM 


T  (pfn  ds 


ASM 

X 


(68) 


where  T{p^^^)  is  given  by  Equation  (49). 

The  ASM  integral  in  Equation  (68)  has  singularities  due  to  the  term  and  the  T-factor. 

Py 

One  can  remove  the  weak  -3W  singularity  with  the  variable  transformation 

Py 


{sf^  +  2k) 

=  sin(i/) 

(69) 

^l-isiSM  +  2kf 

=  cos{u) 

(70) 

dsf^ 

=  cos(i/)  du 

(71) 
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where  +  2k)  =  with  c4  =  ^.  Moving  to  a  shifted  cylindrical  coordinate  system,  one 

defines 


Pt  =  +  y1 

(72) 

(73) 

The  ASM  integral  then  transforms  to 

8in"^(2i:+l) 

A 

(74) 

1) 

where  the  integration  range  traces  out  a  segment  of  the  integration  contour  shown  in  Figure  6.  In 
particular,  for  a  given  k  value,  the  integration  path  lies  on  either  of  the  two  vertical  paths  or  the  real 
axis  between  —  |  and  One  can  handle  the  T-factor  singularities  as  before  by  adding  a  slight  loss 
to  the  denominator  associated  with  propagation  in  the  y  direction. 

Equation  (62),  along  with  Equations  (61),  (68),  and  (74),  gives  the  final  form  for  the  field 
radiated  by  an  edge  element.  As  an  example  of  an  edge  element  calculation,  consider  the  field 
magnitude  plot  in  Figure  19.  As  in  the  PBF  case,  the  edge  element  field  is  sampled  over  an  x  directed 
scan  at  yt  =  O.OSAo  above  the  element.  The  element  is  centered  in  a  0.2Ao  thick  slab  whose  relative 
permittivity  varies  from  1  to  3.25.  The  figure  clearly  shows  the  onslaught  of  surface  waves  as  the 
permittivity  rises,  just  as  Figure  17  did  for  the  PBF. 

2.1.3  Test  Function  Location  Convention.  Sections  2.1.1  and  2.1.2  developed  expressions 
for  the  fields  radiated  by  the  two  basis  function  types.  The  values  of  these  fields  at  the  locations  defined 
by  the  impulsive  testing  functions  are  directly  proportional  to  the  coupling  between  the  particular 
basis  and  testing  functions  through  the  unknown  current  coefficient.  Consequently,  the  author  used 
Equations  (56)  and  (62)  to  calculate  the  coupling  matrix  entries  by  dividing  out  the  unknown  current 
coefficients.  Equation  (9)  gives  the  basis  and  testing  function  ordering  in  the  coupling  matrix.  To 
avoid  confusion.  Figure  20  shows  the  exact  basis  and  testing  function  locations  for  a  semi-infinite 
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Figure  19.  Edge  Element  Radiation  in  the  Presence  of  a  Dielectric  Slab  of  Various  Relative 
Permittivities 

wire  array  with  three  edge  elements.  In  particular,  notice  how  the  PBF  overlaps  the  edge  element 
basis  functions.  In  essence,  the  edge  element  current  coefficients  modify  the  PBF  current  coefficient 
for  the  first  three  elements.  The  author  found  the  optimal  PBF  test  location  lies  one  inter-element 
distance  from  the  inner-most  edge  element  toward  the  infinite  part  of  the  array.  Intuitively,  this  makes 
sense,  since  the  irmer-most  edge  element  should  have  the  strongest  coupling  with  the  infinite  portion 
of  the  array. 

2.2  Excitation  Vector 

Turning  now  to  the  left  hand  side  of  Equation  (9),  consider  the  free-space  geometry  shown  in 
Figure  5.  The  source  of  the  incident  field  lies  out  at  infinity  and  has  an  amplitude  sufficient  to  produce 
a  1.0  (^)  plane  wave  at  the  dielectric  slab’s  left-most  boundary.  Furthermore,  the  Dirac  delta  testing 
functions  convert  the  inner-product  integrals  to  simple  function  evaluations.  The  excitation  vector 
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•  PBF  Element  Location 
°  Edge  Basis  Function  Location 
X  Test  Function  Location 


Figure  20.  Basis  and  Testing  Function  Locations  for  Semi-Infinite  Wire  Array 


entries  are  then 


(75) 


where  p*  is  the  test  function  location  indicated  in  Figure  20. 


One  can  introduce  a  stratified  media  to  the  geometry  using  the  T-factor  concept  once  again. 
Unlike  the  coupling  matrix  case,  the  plane  wave  reference  location  lies  just  outside  the  left-most 
dielectric  boundary.  One  must  first  use  an  effective  transmission  coefficient  to  take  the  plane  wave 
inside  the  slab. 


r, 


12 


_ 712 _ 

1  -f  ri2r23e-j2ft(2'2-3/i)«!,2 


(76) 


where  the  superscript,  e,  identifies  an  effective  Fresnel  coefficient  and  from  Equation  (42) 


7-12  = 


(77) 
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Equation  (76)  takes  into  account  the  reflected  plane  waves  bouncing  between  the  dielectric  boundaries 
to  form  a  composite  field  at  the  left  boundary.  One  obtains  the  T-factor  by  simply  adding  the 
appropriate  phase  to  shift  the  composite  plane  wave  to 

Te,{sy2)  =  T{^e-^My-y.)sy2  (78) 


The  final  expression  for  the  excitation  vector  entries  in  a  dielectric  slab  is  then 

EUp,)  =  (79) 


Note,  one  can  obtain  the  incident  plane  wave  propagation  vector  x  component  by  phase  matching  at 
they  =  yi  boundary. 


5x2  — 


(80) 


2.3  Solution  Vector  and  Far-Field  Radiation 

One  now  has  the  coupling  matrix  and  excitation  vectors  filled.  The  unknown  current  coefficients 
can  be  calculated  from  Equation  (9)  through  a  variety  of  linear  system  techniques.  Plugging  these 
coefficients  into  Equation  (4)  gives  the  approximation  of  the  induced  scattering  currents.  For  the 
particular  basis  function  scheme  shown  in  Figure  20,  one  can  write 

2 

Jz{p)  «  Y.  +  J(6+i))  6  [^{x  -  hd,f  +  y^) 

6=0 

oo 

+ Jo  Y  ^  (V(a:-fe4)2  +  y2)  (81) 

6=3 


This  current  approximation,  plugged  into  the  radiation  integral,  gives  an  expression  for  the  fields 
scattered  by  a  semi-infinite  wire  array.  One  can  then  the  calculate  the  echo-width  from  the  definition 


(72  D  =  lim  27rp 
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(82) 
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Figure  21.  Large  Finite  Wire  Array 


where  \El{p)\  =  1.0  (^)  and  E^{p)  is  calculated  using  the  large  argument  approximation  of  the 
Hankel  function. 

Since,  in  practice,  all  arrays  are  finite,  one  would  like  to  use  the  previous  analysis  to  determine 
the  fields  from  a  large  finite  array.  Consider  the  wire  array  depicted  in  Figure  21.  Using  a  finite 
technique  such  as  Skinner’s [3 3]  to  calculate  the  echo-width,  one  is  faced  with  solving  a  1000  by  1000 
matrix  problem.  Alternatively,  one  could  model  the  large  array  in  two  parts  where  the  currents  on 
the  left  and  right  halves  are  those  computed  for  the  two  corresponding  semi-infinite  arrays.  With  the 
left-most  wire  in  Figure  21  at  the  origin.  Equation  (81)  gives  the  currents  for  the  left  half  elements. 
One  can  calculate  the  right  half  currents  by  solving  the  complement  semi-infinite  wire  array  problem 
shown  in  Figure  22,  where  s  is  reflected  across  the  —  t/  axis.  The  resulting  currents  must  be  given 
the  appropriate  phase  shift  to  place  them  in  the  proper  phase  relationship  with  the  left  half  currents. 
Assuming  three  individual  edge  elements  were  necessary  for  each  semi-infinite  array,  the  combined 
array  problem  results  in  two,  four  by  four  matrix  problems. 
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Figure  22.  Complementary  Semi-Infinite  Array  Problem  for  the  Large  Finite  Array  Right  Edge 

Several  observations  are  now  in  order.  First,  for  most  large  array  problems,  the  system  solution 
eats  up  a  major  portion  of  the  total  computation  time.  As  a  conservative  estimate,  the  number  of 
machine  operations  required  to  solve  the  system  increases  as  with  the  number  of  unknowns.  Thus, 
one  can  achieve  tremendous  computational  savings  using  the  hybrid  PMM/MM  method.  In  addition, 
one  can  increase  the  finite  array  size  in  Figure  21  almost  without  limit  with  virtually  no  computational 
penalty.  This  is  a  consequence  of  the  fact  the  edge  element  perturbations  beyond  the  PBF  test  location 
are  insignificant.  Second,  the  hybrid  PMM/MM  method  has  a  “built-in”  check  for  the  number  of 
required  edge  elements.  Specifically,  if  the  PBF  current  coefficient  calculated  in  the  left  edge  problem 
is  significantly  different  from  the  PBF  current  coefficient  from  the  right  edge  problem,  more  edge 
elements  need  to  be  modeled.  Since  this  requires  a  costly  system  recalculation,  one  should  start  with 
a  good  a-priori  estimate  for  the  number  of  edge  elements.  Hansen  and  Gammon[16]  provide  just 
such  an  estimate  based  on  their  Gibbsian  model  for  finite  scanned  arrays.  Their  estimate. 
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assumes  a  symmetric  edge  element  perturbation  valid  for  incident  angles  away  from  grazing.  Finally, 
note  that  for  the  truly  semi-infinite  array,  the  PBF  does  not  contribute  to  the  far-field  pattern,  except 
perhaps  at  a  finite  number  of  angles  associated  with  grating  lobe  directions.  This  observation  is 
completely  analogous  to  how  the  Physical  Optics  (PO)  currents  do  not  contribute  to  the  far-field 
pattern  in  a  Physical  Theory  of  Diffraction  (PTD)  solution  for  the  scattering  from  a  half-plane. 

This  chapter  developed  the  hybrid  PMM/MM  method  outlined  in  Chapter  I.  Using  a  simple 
two-dimensional  wire  geometry  as  a  proof-of-concept,  the  chapter  clearly  demonstrated  the  method’s 
benefits.  The  actual  validation  of  a  computer  code  based  in  this  development  is  presented  in  Chap¬ 
ter  IV.  The  next  chapter  turns  from  the  proof-of-concept  geometry  to  a  more  practical,  slotted  array 
problem. 
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III.  Three-Dimensional  Semi-Infinite  Slot  Array 


The  sciences  do  not  try  to  explain,  they  hardly  even  try  to  interpret,  they  mainly  make 
models.  By  a  model  is  meant  a  mathematical  construct  which,  with  the  addition  of 
certain  verbal  interpretations,describes  observed  phenomena.  The  justification  of  such  a 
mathematical  construct  is  solely  and  precisely  that  it  is  expected  to  work. 

-  John  von  Neumann 

This  chapter  applies  the  hybrid  PMM/MM  approach  to  the  slotted  ground  plane  shown  in 
Figure  23.  The  array  consists  of  slot  columns  which  can  be  shifted  or  “skewed”  in  the  2  direction. 
A  multi-layer  dielectric  media  encapsulates  the  slotted  array,  providing  a  good  approximation  to  a 
practical  radome  design. 

As  in  the  wire  array  problem,  one  first  examines  the  geometry  with  the  goal  of  reducing  the 
domain  of  the  unknown  currents.  Using  the  equivalence  principle,  one  arrives  at  the  equivalent 
geometry  shown  in  Figure  24.  Note  how  the  now  complete  ground  plane  breaks  the  problem  into 


Figure  23.  Semi-Infinite  Slot  Array  Geometry 
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two  regions.  The  advantage  of  this  particular  equivalent  geometry  lies  in  the  fact  the  only  non-zero 
currents  are  magnetic  currents  restricted  to  the  former  slot  locations.  Equation  (84)  defines  these 
currents  in  terms  of  the  outward  surface  normals  and  the  total  electric  field  in  the  original  geometry. 

-^1,2  =  Exhi^2  (84) 

Even  though  Figure  24  shows  a  thick  ground  plane,  for  this  analysis  the  author  assumes  an  infinitely 
thin  ground  plane.  As  a  result,  the  continuity  of  the  tangential  electric  field  provides  a  coupling 
between  the  two  half-space  problems,  effectively  cutting  the  number  of  unknowns  in  half.  Based 
on  the  tangential  magnetic  field  boundary  condition,  one  can  write  a  magnetic  field  integral  equation 
relating  the  unknown  scattering  currents.  Mi,  to  the  known  incident  fields.  Specifically, 


fiiXiJ"*  =  -fiix  +  Hl{Mi) 


Res 


(85) 
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where  the  generalized  incident  field,  H^\  accounts  for  the  presence  of  the  ground  plane  and  field 
equality  holds  only  at  the  former  slot  locations,  5'.  In  a  manner  analogous  to  the  thin  wire  approx¬ 
imation  of  Chapter  II,  one  can  simplify  Equation  (85)  using  a  thin  slot  approximation.  The  vector 
equation  becomes  a  scalar,  resulting  in  the  relationship  between  the  z  directed  fields, 

Hf  =  -z-\H[(zM,)  +  Hl{zM,)]  (86) 

Equation  (86)  can  be  solved  using  the  PMM/MM  procedure  with  the  current  approximation, 

B 

_  (87) 

6=0  R^es 

and  the  inner  product  definition, 

{9t(R),4>b(R))  =  j  6t{R)4>i[R)dR  (88) 

s 

where  (pbiR)  and  6t{R)  are  periodic  basis  and  testing  functions  associated  with  the  slot  columns. 
Plugging  Equation  (87)  in  Equation  (86)  and  performing  the  inner  product,  one  arrives  at  the  matrix 
equation  analogous  to  Equation  (8)  of  Chapter  I, 

J  Hf(R)et(R)  dR^-J2  Mb  Jz-  [Hl{zcf>b(R!))  +  Hl{zcf>b{R!))\  dR  (89) 

S  6=0  J 

where  t  =  0,1, ...  B.  The  remainder  of  this  chapter  examines  each  component  of  Equation  (89), 
beginning  with  the  coupling  matrix. 

3 . 1  Coupling  Matrix 

The  coupling  matrix  elements  based  on  magnetic  currents  represent  admittances  defined  by 
Equation  (90). 

Ytb  =  -Jz-  ^{{zcl>b(R!))  +  Hl{z(f>b{R!))\  et(R)  dR  (90) 

s 
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Just  as  for  Chapter  IPs  wire  array,  the  basis  functions,  (f)b{R'),  take  on  two  distinct  forms  depending 
on  their  location  in  the  array.  Those  near  the  array’s  edges  model  single  slot  columns,  while  those 
obeying  Floquet  phasing  model  the  entire  center  portion  of  the  array.  The  next  subsection  examines 
the  latter,  begiiming  with  the  homogeneous  media  problem  and  then  introducing  a  stratified  dielectric 
media. 


3.1.1  Physical  Basis  Function. 


3. 1.1.1  Physical  Basis  Function  Radiating  to  Test  Function  in  a  Homogeneous  Media. 

The  slot  array  problem  adds  another  level  of  complexity  to  the  two-dimensional  wire  array  problem 
by  introducing  a  third  dimension.  Fortunately,  the  array  is  fully  periodic  in  this  dimension^,  allowing 
one  to  reduce  the  problem  to  that  of  determining  the  induced  currents  on  a  single  reference  slot.  The 
remaining  currents  in  the  associated  slot  column  are  then  related  to  the  reference  slot  via  a  Floquet 
phase  shift.  In  fact,  looking  at  the  slot  array  as  an  array  of  slot  columns,  one  finds  a  tremendous 
similarity  between  the  wires  of  Chapter  II  and  the  slot  columns  here.  This  fact  allows  one  to  build  on 
the  development  of  Chapter  II  in  deriving  analogous  expressions  for  the  slot  array. 

One  begins  the  analysis  with  the  semi-infinite  array  of  z  directed  Hertzian  dipoles  depicted  in 
Figure  25.  One  can  obtain  the  radiation  from  the  slot  array  by  integrating  the  magnetic  dipole  array 
along  the  slot’s  long  dimension.^  From  Equation  (90)  it  is  apparent  one  needs  the  z  component  of  the 
magnetic  field  radiated  by  the  PBF  in  order  to  calculate  the  admittance.  The  author  uses  the  electric 
vector  potential  to  obtam  the  magnetic  field  produced  by  the  magnetic  currents  [17].  Specifically, 
one  can  write  the  vector  potential  from  the  Hertzian  dipole  as. 


^  e  MpBF  e 


FqmiR)  =  Z 


-j\R-R,m\ 


47r|  J?  —  R, 
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(91) 


where  q  is  the  column  index  associated  with  the  x  direction,  m  is  the  row  index  associated  with 
the  z  direction,  and  Mpsi?  is  the  magnetic  current  value  at  Rqm-  Recalling  the  PBF  obeys  Floquet 


^For  this  research,  the  author  defines  this  to  be  the  z  direction. 
^This  assumes  a  thin  slot  analogous  to  Chapter  II’s  thin  wire. 
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Figure  25.  Hertzian  Dipole  PBF  Radiating  to  Test  Location  in  Homogeneous  Media 


phasing  in  the  semi-infinite  direction,  one  can  express  the  vector  potential  from  the  entire  semi-infinite 
Hertzian  dipole  array  as, 


F{R)  =  z 


€  MpBF  dz‘ 
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E  E 

g=0  m=— 00 
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(92) 


Equation  (92)  represents  a  completely  spatial  domain  expression.  For  reasons  given  in  Chap¬ 
ter  II,  one  desires  a  spectral  domain  form  based  on  a  plane  wave  expansion.  In  the  z  direction,  one 
can  use  the  Poisson  sum  formula  to  transform  the  row  summation  to  the  spectral  domain.  Plugging 
the  transform  pair  [3], 


+  {(jj  —  zy 


(93) 
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in  Equation  (174)  of  Appendix  A,  one  obtains, 


F(R) 


^  €  Mpgp  dz' 
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E 


^=0 


n^z\  \ 
dx  dz  ' 


(94) 


where, 


(95) 


Equation  (94)  is  a  mixed  spectral/ spatial  domain  form,  hereafter  referred  to  as  the  “spatial” 
domain  expression.  Notice  the  similarities  between  the  inner  summation  in  Equation  (94)  and  the 
spatial  domain  expression  for  the  wire  in  Equation  (17).  Clearly,  the  slot  column  radiates  cylindrical 
waves  which  propagate  or  become  evanescent  depending  on  the  direction  component,  r^. 

To  obtain  the  complete  spectral  domain  expression,  one  uses  the  one-sided  Poisson  sum  formula 
given  in  Equation  (185)  to  transform  the  column  summation  in  Equation  (94)  to  the  spectral  domain. 
With  the  Fourier  transform  pair  given  by  Equation  (18)  and  a  few  algebraic  manipulations,  one 
obtains. 
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where, 


nAzA  k\ 

^  dxdz  dx 


1 7*2  —  /j-2 

p  ’  X 


As  for  the  wire  array,  Equation  (96)  is  a  plane  wave  expansion  containing  both  discrete  and  continuous 
spectra. 

Equations  (92),  (94),  and  (96)  give  three  alternate  forms  of  the  electric  vector  potential  radiated 
by  the  Hertzian  dipole  array.  Using  these  expressions  along  with  the  superposition  principle  and 
image  theory,  one  can  obtain  the  vector  potential  from  a  finite  length  magnetic  dipole  array  in  the 
presence  of  a  ground  plane.  The  thin  slot  assumption,  applied  to  z  directed  slots,  converts  the 
superposition  integral  into  a  line  integral  in  z.  If  one  assumes  the  basis  function  has  a  piecewise 
sinusoidal  (PWS)  shape. 


Sin  (^) 


(100) 


one  can  express  the  spectral  domain  form  of  the  vector  potential  from  either  region  1  or  region  2  as. 
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where  zi,  is  the  center  of  the  PWS  and  is  the  integral  evaluation, 


pb  ^  2[cos(^)  -  cos(^)]  ^ 
/?sin  (f)il-rl) 


referred  to  as  a  “pattern  factor.”  Equiped  with  the  vector  potential,  one  can  now  determine  the 
magnetic  field  from  the  relationships  [17], 

E\R)  =  ^VxF(R)  (103) 

H\R)  =  —VxWCR)  (104) 

JUfi, 

Plugging  Equation  (101)  into  Equation  (103)  gives  the  spectral  domain  form  of  the  scattered  electric 
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and  Equation  (105)  into  Equation  (104)  gives  the  z  component  of  the  scattered  magnetic  field. 
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where  V  is  the  intrinsic  admittance  of  the  homogeneous  media.  The  spatial  domain  forms  are  found 
by  plugging  Equation  (94)  into  Equation  (103), 
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and  Equation  (107)  into  Equation  (104), 
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Equations  (106)  and  (108)  give  the  z  component  of  the  magnetic  field  radiated  by  the  PBF  in  a 
homogeneous  media.  One  obtains  the  coupling  matrix  terms  associated  with  the  PBF  by  taking  the 
inner  product  of  either  Equation  (106)  or  Equation  (108)  with  a  PWS  testing  function,  as 

indicated  by  Equation  (90).  Before  proceeding,  however,  it  is  numerically  expedient  to  transform 
Equation  (106)  to  the  angular  spectral  domain. 

Defining  the  continuous  spectral  integral  as, 

7  f  +  y^/rl  -  tI] 

"  /  v<-'i  -  -  ’■.) 

—  OO 


one  must  consider  two  cases  depending  on  the  value  of  Vp.  The  first  case  occurs  when  <  1,  making 
Tp  a  strictly  positive,  real  number.  One  can  then  divide  by  r^,  resulting  in. 
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(110) 


where  r'  =  ^.  Equation  (110)  now  has  exactly  the  same  form  as  the  spectral  integral  in  Equa¬ 
tion  (21).  As  before,  one  removes  the  weak  singularity,  ,  ^  from  the  integrand  via  the  variable 
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substitution,  =  sin  {u).  With  the  displaced  coordinate  system,  p  =  [x  +  and 

(f)  =  tan“^  I  ^  ^  '  I ,  one  obtains  the  angular  spectral  integral. 


4  (sin(r.)-^) 


-jPr^p  cos  V 


(111) 


where  the  integration  contour,  C,  is  shown  in  Figure  6.  Finally,  one  extracts  the  singularity,  Uo  = 
sin“^  (?^)  ’  resulting  in  the  expression. 
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and. 


where  one  chooses  either  the  upper  or  lower  line  for  j  <  1  or  j  >  1,  respectively. 

The  second  case  occurs  when  >  1.  Using  the  variable  substitution,  fp  =  y/rl  —  1,  the 
weak  singularity  goes  away,  removing  the  need  to  transform  to  the  angular  spectrum.  The  continuous 
spectral  integral  can  now  be  written. 


Extracting  the  singularity,  one  has. 


where. 


(120) 
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Note  the  singular  integrand  is  an  odd  function  causing  the  integral  to  evaluate  to  zero.  Equations  (112) 
and  (1 19)  are  the  continuous  spectral  forms  the  author  implemented  in  the  computer  codes  validated 
in  Chapter  IV. 

At  this  point,  one  may  proceed  with  the  inner  product  integral  to  obtain  the  desired  admittance 
matrix  element.  From  Equation  (90),  there  are  actually  two  integrations  to  perform.  These  “one¬ 
sided”  admittances  correspond  to  the  two  half-spaces  created  by  the  ground  plane.  For  symmetric 
dielectric  stacks,  the  admittances  are  identical,  requiring  one  to  compute  only  one  side  and  double  the 
results.  The  spatial  domain  form  comes  from  Equation  (108).  Assuming  a  PWS  test  function  defined 
to  be  non-zero  over  the  reference  slot  centered  atR  =  Rt,  one  can  write. 
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where  the  basis  and  testing  function  pattern  factors  are. 
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In  a  similar  manner,  the  spectral  domain  form  comes  from  Equation  (106), 
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Figure  26.  Geometry  for  PBF  One-sided  Admittance  Problem 

In  order  to  avoid  singularities  in  the  self  admittance  calculations,  the  basis  and  testing  functions  must 
not  coincide.  As  before,  one  places  the  testing  function  an  “effective  wire  radius”  out  of  the  array 
plane.  For  thin  slots,  this  is  defined  to  be  a  =  slot  .^idth 

This  section  developed  both  the  spatial  and  spectral  domain  forms  for  the  coupling  between 
a  finite  length  periodic  basis  function  and  a  single  finite  length  testing  function.  The  homogeneous 
problem  again  forms  the  foundation  for  extending  the  development  to  a  stratified  dielectric  media. 
The  next  section  uses  Equations  (121)  and  (123)  to  develop  the  admittances  for  the  slot  array  encased 
in  a  stratified  dielectric  media. 

3. 1.1. 2  Physical  Basis  Function  Radiating  to  Test  Function  in  a  Stratified  Media. 
Figure  26  depicts  the  one-sided  admittance  problem  involving  the  PBF  radiating  to  a  test  location  in 
the  presence  of  a  stratified  dielectric  media.  Similar  to  the  wire  array,  one  can  use  the  T-factor  concept 
to  account  for  the  presence  of  the  dielectric  layers.  The  plane  wave  nature  of  the  spectral  domain 
expressions  allows  one  to  break  the  propagation  into  a  direct  term  and  a  term  accounting  for  the 
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multiple  bounces  between  the  dielectric  layers.  Unlike  the  wire  array,  Figure  26  shows  a  multi-layer 
geometry.  Consequently,  the  T-factor  expression  used  here  is  slightly  more  complicated. 

Recall  the  author  uses  Munk’s  convention  to  define  the  effective  reflection  and  transmission 
coefficients.  Also  note  the  expressions  developed  heretofore  assume  the  test  locations  have  a  positive 
y  component^.  Consequently,  the  z  polarized  magnetic  fields  require  the  following  magnetic  field 
related  coefficients. 


t'M+1  = 
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'n(i+i)ry{i)  -b  V(i)'>'y(i+i) 
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where  the  e  superscript  denotes  an  effective  coefficient.  Equations  (124)  and  (125)  are  Fresnel 
coefficients  which  give  the  transmitted  and  reflected  field  ratios  at  the  dielectric  interface  between 
two  half-spaces.  The  effective  coefficients  give  the  composite  field  ratios,  accounting  for  the  multiple 
reflections  and  transmissions  in  the  multi-layer  structure.  The  coefficient  for  a  particular  layer  depends 
on  those  from  the  next  layer  in  the  y  direction.  Therefore,  one  must  start  with  the  outer  half-space 
layer  and  iteratively  compute  the  coefficients  down  to  the  layer  closest  to  the  ground  plane.  Equiped 
with  this  inner-most  effective  reflection  coefficient,  one  can  repeat  the  T-factor  analysis  performed  in 
Chapter  II. 

^This  restriction  does  not  limit  the  analysis.  To  account  for  region  1,  one  simply  flips  the  geometry  and  mirrors  the 
excitation  vector  through  the  z  axis. 
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In  particular,  one  can  break  the  plane  waves  into  four  summations  as  shown  in  Figure  13.  These 
can  be  written  in  closed  form  and  combined  to  give  an  expression  for  the  field  at  Rt, 


^direct 
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where  is  the  outer  boundary  of  the  layer,  yj  =  0,  and  Tfi_i  =  1.  Expression  (128)  has  a  direct 
and  bounce  component.  Note  however,  the  homogeneous  expressions  developed  in  the  previous 
section  actually  double  the  direct  term  through  image  theory  in  order  to  account  for  the  ground  plane. 
To  use  Expression  (128),  one  must  pull  out  a  factor  of  two  resulting  in, 
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where  the  second  term  in  the  square  brackets  is  defined  to  be  the  T-factor, 
Equiped  with  the  T-factor,  one  may  now  rewrite  Equation  (106), 
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The  pattern  factor  here  is  modified  in  accordance  with  Kombau’s  [21]  work  to  account  for  the 
possibility  of  different  permittivities  on  either  side  of  the  array  ground  plane.  Kombau  used  the 
arithmetic  average  of  the  permittivities  in  his  calculations.  Here,  the  author  found  the  arithmetic 
average  of  the  wave  numbers  worked  equally  as  well.  The  redefined  pattern  factor  is  then. 
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where  Pavg  =  The  continuous  spectral  integral  can  be  expressed. 
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where, 
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and  the  analytic  solutions  for  the  extracted  integrals  are  given  in  Equations  (116)  and  (117),  or 
Case  >  1: 
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The  previous  development  gives  the  fields  radiated  to  the  test  location  minus  the  direct  term. 
To  arrive  at  the  final  expression  for  the  admittance  between  the  PBF  and  the  test  function,  one  simply 
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combines  the  spatial  domain  expression  in  Equation  (121)  for  the  direct  term  with  Equation  (130)  for 
the  bounce  term  to  obtain, 
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where  P/  is  the  test  function  pattern  factor  obtained  by  changing  the  sign  on  the  exponent  in 
Equation  (131).  The  author  implemented  Equation  (138)  in  the  slot  code  validated  in  Chapter  IV.  The 
T-factor  singularities  are  handled  as  before  by  introducing  a  slight  dielectric  loss  in  the  y  propagation 
direction  and  increasing  the  integrand  sampling  in  the  neighborhood  of  the  former  singularity. 

This  completes  the  development  of  the  coupling  matrix  terms  associated  with  the  physical  basis 
function.  The  next  section  derives  expressions  for  the  coupling  terms  associated  with  the  individual 
edge  element  basis  functions.  With  these,  the  coupling  matrix  is  complete. 


3.1.2  Edge  Element  Basis  Function. 


3. 1.2.1  Edge  Element  Basis  Function  Radiating  to  Test  Function  in  a  Homogeneous 
Media.  Consider  a  single  periodic  basis  function  entered  at  R},.  One  can  obtain  the  coupling 
between  this  column  and  a  reference  test  mode  located  at  Rt  from  Equation  (121), 
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where  P^  and  P*  are  defined  by  Equation  (122).  Recall  this  spatial  domain  expression  results  from 
the  application  of  image  theory  to  replace  the  array  ground  plane  with  a  homogeneous  media.  Adding 
a  stratified  dielectric  media  complicates  the  problem  considerably.  Once  again,  the  author  chose  to 
take  advantage  of  the  plane  wave  expansion  in  the  ASM  to  handle  the  stratified  media. 
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3. 1.2. 2  Edge  Element  Basis  Function  Radiating  to  Test  Function  in  a  Stratified  Media. 

The  basic  ASM  approach  involves  creating  a  pseudo-array  around  the  edge  element  in  order  to 
transform  to  the  spectral  domain.  One  can  then  use  the  T-factor  developed  in  the  PBF  section  to 
account  for  the  stratified  media.  One  obtains  the  fields  from  the  original  edge  element  via  a  Fourier 
transform.  Defining  the  pseudo-array  inner-column  spacing,  ,  such  that  -1-  2k, 

one  obtains  the  following  expression  for  the  pseudo-array  coupling, 
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where  write  the  coupling  from  the  single  edge  element  in  terms 

of  the  ASM  integral  as. 
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where  the  bounce  superscript  indicates  the  direct  coupling  is  not  included.  Next,  one  finds  it 
numerically  expedient  to  transform  to  the  angular  spectral  domain. 
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where. 
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All  that  remains  is  to  combine  the  bounce  term  in  Equation  (142)  with  the  direct  term  in  Equation  (139)'^ 
to  arrive  at  the  coupling  expression  implemented  in  the  author’s  code. 

Ytb  =  (145) 

Equations  (138)  and  (145)  define  all  the  entries  in  the  coupling  matrix.  In  order  to  expedite  the 
edge  treatment  study,  the  author  also  incorporated  the  ability  to  apply  a  complex  impedance  to  the 
individual  slot  columns.  The  next  section  details  the  required  analysis. 

3.1.3  Edge  Element  Loading.  Consider  the  reference  element  of  a  periodic  basis  function. 
A  typical  slot  load  consists  of  a  distributed  resistance  such  an  “r-card.”  One  would  like  to  convert 
the  distributed  r-card  resistance  to  a  lumped  admittance  which  could  then  be  added  to  the  diagonal 
elements  of  the  coupling  matrix.  Figure  27  depicts  this  concept  for  the  PWS  magnetic  current  mode. 
The  author  used  a  method  suggested  by  English  [9]  to  perform  the  conversion.  One  simply  equates 
the  dissipated  power  in  both  representations.  From  Equation  (84),  an  electric  field  perpendicular  to 
the  slot’s  length  produces  a  z  directed  magnetic  current.  Thus,  one  can  express  the  power  in  the 
distributed  representation  as, 

Po  =  J  dP 

slot 


L 

2 


2 

L 

^  RpM^  }  [sin  {l3{^-\z\)) 

2LArJ^  [  sin(f) 

^  All  quantities  associated  with  the  material  parameters  are  taken  to  be  those  of  the  inner-most  layer,  1. 
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where  AR  is  the  slot’s  aspect  ratio.  The  power  in  the  lumped  representation  is, 


Po  = 


RlMI 

2 


(147) 


Equating  the  two  and  solving  for  Rl,  one  obtains. 


AR 


PL-  sin  ipL) 
2/3Lsin  (^) 


(148) 


The  reciprocal  of  Equation  (148)  gives  the  lumped  admittance  which  the  author  adds  to  the  appropriate 
self-admittance  term  in  the  coupling  matrix. 

The  previous  development  assumes  the  mode  represents  the  total  edge  element  current  rather 
than  a  perturbation  of  the  PBF  current,  as  in  the  scheme  depicted  in  Figure  20.  It  is  possible  to  derive 
the  total  edge  element  formulation  from  the  perturbation  formulation  for  an  unloaded  array.  One 
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PBF  Test  Location 
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I 

•  PBF  Element  Location 
°  Edge  Basis  Function  Location 
X  Test  Function  Location 

Figure  28.  Alternate  Basis  and  Testing  Function  Locations 

can  then  apply  the  edge  element  loading  to  the  new  coupling  matrix.  Skinner  suggests  the  following 
method  to  convert  between  the  two  unloaded  formulations  [32]. 

First,  consider  the  perturbation  formulation’s  matrix  equation, 


Tl,!  ^1,2  • 

T2,1  ••• 

yi,N  Yi,N+1 

y2,N  y2,N+l 

^  '-5^  ■' 

1 _ 

Ym,i  •  •  • 

Yn^n  kjv,iv+l 

Vn 

In 

_  Liv+l,!  •  •  • 

Tiv+i,iv  Yjv_|_i,iv’+i 

Yn+i 

In+1 

where  the  subscripts  Ito  N  denote  edge  elements  and  N  +  1  denotes  the  PBF.  The  PBF  test  location 
is  indicated  in  Figure  20. 

Now,  consider  sliding  the  PBF  to  the  right  so  it’s  left  edge  element  coincides  with  the  old 
PBF  test  location,  as  in  Figure  28.  The  edge  modes  in  the  new  formulation  represent  the  entire  edge 
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element  currents.  Keeping  the  old  PBF  test  location,  one  can  express  the  new  matrix  equation  as 


Yi,i  n,2  •••  h 

^2,1  y2,N  y2,N+l  ^2  h 

\  \  \  \  =  \  (150) 

yN,i  •  ■  ■  yN,N  yN,N+i  Vat 

^iV+1,1  •  •  •  yN+l,N  ^JV+l,iV+l  l^r+l  ^JV+1 

where  the  tilde  indicates  a  different  value  than  the  perturbation  formulation. 

The  key  to  Skiimer’s  method  lies  in  recognizing  the  following  equivalences  between  the  two 
formulations,  based  on  the  assumption  both  yield  the  same  magnetic  currents. 


Vi  +  Vjv+i 
V2  +  VV+i 

Vn  +  Vat+i 
Vn+1 


(151) 


Once  one  has  solved  the  perturbed  system,  the  unknown  unperturbed  coupling  matrix  entries  can 
be  found  using  the  new  currents  from  Equation  (151)  and  the  known  coupling  matrix  entries  from 
Equation  (149).  For  example,  solving  the  first  row  in  Equation  (150)  for  the  unknown  coupling  matrix 
entry  yi,jv+i  yields  the  equation. 


h  —  yi,iVi  —  -  yi,NyN 


where  the  right  hand  side  values  are  known.  Equiped  with  the  new  coupling  matrix  representing  total 
currents,  one  can  now  apply  the  shunt  admittances  to  the  diagonal  entries  as  described  earlier. 
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Figure  29.  Plane  Wave  Illumination  of  Slot  Column 


3.2  Excitation  Vector 

The  left  hand  side  of  Equation  (89)  defines  the  excitation  vector  entries, 

J  Hf(R)9t(R)  dR  (153) 

s 

where  the  generalized  incident  field  accounts  for  the  presence  of  the  ground  plane.  Since  a  plane 
wave  illuminates  the  structure,  the  author  once  again  uses  a  T-factor  to  account  for  the  dielectric 
layers  shown  in  Figure  29.  One  has  complete  freedom  to  define  the  incident  plane  wave.  Since  the 
validation  data  is  entirely  monostatic,  the  author  chose  to  restrict  the  incident  plane  wave  propagation 
to  the  xy  plane.  In  particular,  the  author  chose  a  l(£),  z  polarized  magnetic  field  with  a  phase 
reference  at  Rg.  Assuming  a  homogeneous  media,  the  magnetic  field  at  the  test  location,  (a;*,  0,  Zt), 
is  then, 

{^R^  —  £~i9o(xiSxo  "F  dfQfdiSyo)  (154) 
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Introducing  the  stratified  media,  one  can  find  the  y  propagating  field  just  inside  the  first  slab  from 
Equation  (126), 


_ _ 


(155) 


One  then  propagates  the  fields  through  the  remaining  layers  to  the  test  location  Rt. 


H'ARt) 


7-1 


U=i 


(156) 


To  get  the  total  field  at  Rt,  one  sums  the  composite  right-going  field. 


Bi(Rt) 


7-1 


.(=1 


and  the  composite  left-going  field. 


(157) 
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(158) 


to  arrive  at. 


HiiRt) 


■  /-I 

2  JJ  r|_i 
.  ^=1 


Hi{Rt)Tex{Sxo) 


(159) 


where  Ti_ii  =  1.  To  account  for  the  irmer  product  with  the  testing  function,  one  uses  the  same 
pattern  factor  as  before,  arriving  at. 


J  HI^{R)et{R)  dR  =  T,,{s,o)  PU  (160) 

5 


3.3  Solution  Vector  and  Far-Field  Radiation 

With  the  coupling  matrix  and  excitation  vectors  filled,  one  can  solve  the  system  to  obtain 
the  unknown  expansion  coefficients.  One  simply  plugs  these  coefficients  back  into  the  current 
approximation  to  obtain  the  desired  scattering  currents.  For  the  code  validated  in  Chapter  IV,  the 
author  defined  the  edge  element  currents  as  perturbations  of  the  PBF  current.  Thus,  the  scattering 
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current  amplitudes,  Mb,  discussed  below  are  actually  a  sum  of  the  PBF  and  the  particular  edge 
element. 

First,  consider  the  homogeneous  media  case  involving  a  periodic  basis  function  centered  at 
{xb,  0,  Zb).  The  associated  scattered  magnetic  field®  can  be  written. 


Ht{R) 


—ojeiMb 

2d^ 


E 


Vi 


(161) 


where  pb  =  yj {xb  —  x)^  +  y'^ .  Restricting  the  plane  of  incidence  to  the  xy  plane,  one  is  interested 
in  the  scattered  field  a  large  distance  away  from  the  array.  For  monostatic  scattering,  only  the  n  =  0 
term  contributes  to  the  far-field.  One  can  use  the  large  argument  form  of  the  Hankel  function, 

to  obtain  a  far-field  approximation  of  the  scattered  magnetic  field, 

~  ™ M)pi  (163) 

2d,  Y  Tt^iTpip 

where  p  =  y/ x'^  -\-y^  and  r]oba  is  measured  from  the  negative  y  axis  in  the  plane  of  incidence.  Note, 
Equation  (163)  contains  the  approximations  pt  ~  p  and  pt  «  p  +  x^sin  (po6«)  for  the  amplitude 
and  phase,  respectively.  One  then  sums  the  field  contributions  from  each  periodic  basis  function®  to 
obtain  the  total  scattered  field, 

®One  only  needs  the  z  component  since  all  the  reference  data  is  co-polarized. 

®  For  the  central  portion  of  a  large  array,  one  simply  adds  the  appropriate  phase  shift  to  the  PBF  coefficient  and  applies 
the  result  to  the  corresponding  slot  column. 
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In  order  to  account  for  the  stratified  media,  the  author  employed  the  ASM,  using  the  T-factor 
defined  by  Equation  (159).  One  can  then  write  the  scattered  field  from  a  single  periodic  basis  function. 


Ht{R)  = 


-Mb 

riid.Xi 


E  E  pt  (165) 

n——aok——oo  - 


Once  again  keeping  only  the  central  term  and  restricting  the  plane  of  incidence  to  the  xy  plane,  one 


obtains, 


£/•«  rpx  I  c _  (  ASM)  J^ASM 


In  order  to  remove  the  weak  singularity,  one  finds  it  desirable  to  transform  to  the  angular  spectral 
domain. 


T]id^Xi  j 


Tex{s\n{v))  du 


where. 


Pb  =  Y  (a:  -  Xbf  +  2/2 
^b  =  tan-i(2f^) 


(168) 

(169) 


The  integrand  in  Equation  (167)  is  highly  oscillatory  for  large  pb  values,  making  numerical 
integration  virtually  impossible.  Fortunately,  the  integral  lends  itself  to  a  stationary  phase  evaluation. 
In  particular,  one  can  identify  (dipb  as  the  large  parameter  and  —  4>b  as  the  stationary  phase  point. 
With  the  resulting  asymptotic  integral  evaluation,  one  can  express  the  scattered  field  as. 


(170) 


Summing  over  the  modes  and  applying  the  same  approximations  to  pb  as  in  the  homogeneous  case, 
one  arrives  at, 


_  —e^Pip 

KffiR) 
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(171) 


p— +00 


where  sin  (i7o)  is  approximated  by  ^  sin  (77064)* 

In  order  to  compare  with  the  measured  data,  the  author  uses  the  echo  width  definition. 


HiiR)? 

(72d  =  lim  27rp - 

Hi{R)\ 


and  Equation  (171)  to  arrive  at  the  expression, 

B 

^  Mb  sin  (77064)  T,,  sin  (77064))  Pf 


(172) 
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(173) 


The  author  implemented  Equation  (173)  in  the  computer  code  validated  in  Chapter  IV.  In  particular, 
the  author  used  the  left  and  right  semi-infinite  array  approach  described  in  Chapter  II  to  model  the 
large  finite  slot  arrays  in  the  next  chapter. 

This  chapter  applied  the  hybrid  PMM/MM  method  to  a  slotted  ground  plane  embedded  in 
a  general  stratified  dielectric  media.  By  modeling  the  individual  slot  colunms  with  periodic  basis 
functions,  the  author  was  able  to  build  on  the  wire  array  analysis.  The  next  chapter  examines  the 
computer  codes  based  on  this  development  in  order  to  validate  the  method. 
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IV.  Validation 

The  triumphant  vindication  of  bold  theories  -  are  these  not  the  pride  and  justification  of 

our  life’s  work? 

-  Sherlock  Holmes  (Sir  Arthur  Conan  Doyle) 

The  previous  two  chapters  applied  the  hybrid  PMM/MM  technique  to  a  two-dimensional  wire 
array  and  a  three-dimensional  slot  array.  This  chapter  validates  the  hybrid  technique  using  a  two 
pronged  approach  based  on  both  reference  code  and  measurement  comparisons.  The  reference  codes 
allow  one  to  directly  examine  the  currents  induced  on  the  array  elements.  As  mentioned  earlier,  this 
approach  provides  a  more  rigorous  test  of  the  technique’s  accuracy  by  avoiding  the  radiation  integral’s 
smoothing  effect.  Measurement  comparisons  are  included  despite  the  complications  associated  with 
simulating  an  infinite  ground  plane  because,  in  many  respects,  “the  proof  is  in  the  pudding.”  The 
author  felt  the  best  way  to  demonstrate  the  hybrid  PMM/MM  technique’s  usefulness  is  to  test  its 
ability  to  model  “real  world”  radomes.  As  a  consequence  of  the  validation  approach,  the  chapter 
breaks  into  two  parts  corresponding  to  the  reference  code  and  measurement  comparisons.  Within 
each  part,  the  subject  wire  and  slot  array  geometries  are  examined  with  an  eye  to  not  only  build 
confidence  in  the  technique,  but  also  to  highlight  the  situations  where  it  is  most  useful.  This  chapter 
concludes  with  a  motivating  illustration  of  the  hybrid  technique’s  potential  for  handling  large  arrays. 

4.1  Reference  Code  Comparison 

One  has  a  distinct  advantage  when  using  a  reference  code  to  validate  a  new  technique.  Unlike 
measurements,  one  has  tremendous  flexibility  in  choosing  representative  array  geometries  to  truly 
“wring  out”  the  technique.  For  the  hybrid  PMM/MM  technique,  the  author  attempted  to  find  array 
geometries  large  enough  to  ensure  the  central  element  currents  obey  Floquet  theory,  yet  small  enough 
to  allow  the  MM  based  reference  codes  to  solve  the  problems  in  a  reasonable  time. 

In  addition  to  choosing  appropriate  problem  geometries,  one  must  be  careful  in  choosing 
and  using  the  reference  codes.  The  reference  code  provides  the  “truth”  baseline  against  which  one 
measures  the  new  technique.  Obviously,  one  must  have  a  certain  amount  of  faith  in  the  reference 
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(600  Elements) 


Figure  30.  600  Element  Wire  Array  in  Homogeneous  Media  Illuminated  by  Plane  Wave  at  an 
Arbitrary  Incident  Angle 


code’s  accuracy.  The  author  made  every  attempt  to  use  independently  validated  reference  codes.  In 
the  following  material,  cases  where  the  reference  code  results  may  be  questionable  are  clearly  noted. 


4.1.1  Wire  Array  in  a  Homogeneous  Media.  The  two-dimensional  wire  array  serves  as  the 
prototype  geometry  throughout  this  research.  Figure  30  depicts  the  specific  scattering  scenario  the 
author  chose  to  validate  the  hybrid  PMM/MM  technique.  This  particular  array  geometry  represents 
a  typical  example  of  many  tests  the  author  performed  in  developing  the  hybrid  technique. 

The  next  series  of  figures  compares  the  induced  scattering  current  magnitude  and  phase  for  each 
element  in  the  600  element  array.  The  reference  values  are  generated  by  a  standard  600  unknown  MM 
analysis  based  on  the  two-dimensional  free-space  Green’s  function.  For  clarity,  the  figures  leave  out 
the  central  portion  of  the  array.  Enough  elements  are  included  however,  to  show  where  the  currents 
are  behaving  according  to  Floquet  theory. 
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Figures  31  and  32  show  the  induced  currents  on  the  left  edge.  Although  not  shown,  the  right 
array  edge  is  a  mirror  copy  of  the  left  edge.  The  reference  values  are  denoted  by  Xs  and  the  hybrid 
PMM/MM  values  are  denoted  by  Os.  The  array  is  illuminated  by  a  plane  wave  incident  normal  to 
the  array  plane.  The  hybrid  technique  used  20  edge  elements  for  both  left  and  right  hand  problems. 
These  figures  clearly  demonstrate  the  validity  of  the  hybrid  technique. 

Figures  33  through  36  show  equally  as  valid  results  for  the  array  illuminated  5°  off  normal. 
In  this  scattering  problem,  the  author  used  30  edge  elements  for  the  leading  edge  (left  edge)  and 
20  elements  for  the  trailing  edge  (right  edge).  Note  the  linear  phasing  from  element-to-element 
indicating  the  onset  of  the  Floquet  region. 

Figures  37  through  40  show  the  array  illuminated  45°  off  normal.  Clearly,  the  leading  edge 
perturbations  are  more  pronounced.  The  author  again  used  30  edge  elements  for  the  left  edge  hybrid 
PMM/MM  problem  and  20  edge  elements  for  the  right.  Note  the  larger  element-to-element  phase 
rotation  caused  by  the  off  normal  illumination.  Once  again,  the  hybrid  technique  captures  the  correct 
induced  current  behavior.  In  particular,  note  the  smooth  transition  from  the  inner-most  edge  element 
to  the  PBF  elements.  The  PBF  coupling  mitigates  the  edge  effects  normally  present  in  a  finite  array 
made  up  of  the  edge  elements,  providing  a  distinct  improvement  over  the  hybrid  approach  proposed 
by  Cwik  and  Mittra  [8]. 

Figures  41  through  44  show  the  currents  induced  by  a  near  grazing  illumination  of  85°  off 
normal.  The  author  used  100  edge  elements  for  the  leading  edge  and  3  edge  elements  for  the  trailing 
edge  in  this  problem.  Increasing  the  number  of  trailing  edge  elements  does  little  to  improve  the 
solution  accuracy.  The  author  traced  the  problem  to  the  direct  coupling  term  in  the  PBF.  Specifically, 
convergence  inaccuracies  in  the  spatial  domain  summation  acceleration  appear  to  be  the  culprit.  For 
this  extreme  case,  the  author’s  implementation  of  the  hybrid  PMM/MM  technique  begins  to  have 
difficulty  for  the  trailing  edge  element  behavior.  This  trend  also  appears  in  other  plots  for  angles  very 
near  grazing.  As  discussed  in  the  conclusion,  a  different  acceleration  algorithm  should  eliminate  this 
difficulty. 
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Incident  Angle 

MM  Code 
(Unknowns) 

MM  Code 
(cpu  secs) 

PMM/MM  Code 
(Unknowns,  Left/Right  Problem) 

PMM/MM  Code 
(cpu  secs) 

0° 

35.6 

0.5 

5° 

35.2 

30/20 

0.6 

45° 

600 

35.2 

30/20 

1.2 

OO 

o 

600 

35.5 

2.4 

Table  1 .  Wire  Array  in  Homogeneous  Media  Reference  MM  Code  versus  Hybrid  PMM/MM  Code 
Time  Comparison 

To  conclude  this  subsection  on  the  wire  array  in  a  homogeneous  media,  Table  1  highlights  one 
of  the  hybrid  PMM/MM  technique’s  benefits.  The  author  measured  the  execution  times  required  for 
the  MM  and  hybrid  PMM/MM  computer  codes  to  calculate  each  of  the  scattering  problems  on  a  Sun 
Microsystems  Sparc20  Workstation.  Although  not  strictly  an  “apples-to-apples”  comparison  due  to 
the  brute  force  integration  algorithms  the  author  used  in  the  hybrid  code,  the  CPU  times  give  a  rough 
estimate  of  the  computational  savings  possible.  This  information,  coupled  with  the  fact  the  MM 
computational  burden  scales  with  system  size  as  where  the  hybrid  PMM/MM  approach  scales  as 
clearly  demonstrates  the  PBF  approach’s  value  in  large  problems^.  Note  that  unlike  the  MM,  the 
hybrid  technique  does  require  recalculation  of  one  row  in  the  coupling  matrix  for  each  monostatic 
angle.  Thus,  in  generating  monostatic  patterns,  the  hybrid  method  scales  as  and  may  not  provide 
as  dramatic  results,  depending  on  the  problem  size.  Of  course,  the  larger  the  problem,  the  better 
the  hybrid  technique’s  performance  relative  to  the  MM.  Also,  no  coupling  matrix  recalculation  is 
required  for  bistatic  patterns. 


^This  is  true  as  long  as  the  perturbations  of  the  array’s  edge  elements  are  adequately  modeled  by  the  number  of  edge 
elements  used  in  the  hybrid  PMM/MM  analysis. 
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Figure  3 1 .  Induced  Current  Magnitude  on  Left  50  Elements  ofWire  Array  with  Normal  Illumination 


Element  Number 


Figure  32.  Induced  Current  Phase  on  Left  50  Elements  of  Wire  Array  with  Normal  Illumination 
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Element  Number 
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Figure  35.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  with  5°  off  Normal 
Illumination 
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Figure  36.  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  with  5°  off  Normal 
Illumination 
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Element  Number 


Figure  39.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  with  45°  off  Normal 
Illumination 


Figure  40.  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  with  45°  off  Normal 
Illumination 
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o  Hybrid  PMM/MM  Solution 
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Figure  41.  Induced  Current  Magnitude  on  Left  150  Elements  of  Wire  Array  with  85°  off  Normal 
Illumination 
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Figure  42.  Induced  Current  Phase  on  Left  150  Elements  of  Wire  Array  with  85°  off  Normal 
Illumination 
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Figure  43.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  with  85°  off  Normal 
Illumination 
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Figure  44.  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  with  85°  off  Normal 
Illumination 
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Figure  45.  600  Element  Wire  Array  in  a  Dielectric  Slab  Illuminated  by  Plane  Wave  at  an  Arbitrary 
Incident  Angle 

4.1.2  Wire  Array  in  a  Stratified  Media.  This  subsection  turns  to  the  problem  of  a  wire 
array  embedded  in  the  dielectric  slab  shown  in  Figure  45.  The  author  chose  a  MM  code  based  on  the 
ASM  to  provide  the  reference  solutions.  Again,  the  wire  array  contains  600  elements  of  which  only 
the  left  and  right  edges  are  displayed. 

Figures  46  and  47  show  the  wire  array’s  left  edge  subject  to  normal  incident  illumination.  As 
in  the  homogeneous  case,  the  induced  currents  are  symmetric  about  the  middle  of  the  array.  The 
author  used  only  20  edge  elements  as  before.  Here  the  slab  has  strengthened  the  coupling  between 
the  elements,  necessitating  the  inclusion  of  more  edge  elements  to  completely  capture  the  edge 
perturbations.  Note  however,  the  elements  that  were  modeled  are  very  accurate. 

Figures  48  through  5 1  show  the  results  for  the  array  illuminated  5°  off  normal.  Here  the  author 
used  50  and  30  edge  elements  respectively  in  the  hybrid  PMM/MM  analysis.  Once  again,  the  hybrid 
technique  accurately  models  the  induced  currents. 
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Incident  Angle 

MM  Code 
(Unknowns) 

MM  Code 
(cpu  secs) 

PMM/MM  Code 
(Unknowns,  Left/Right  Problem) 

PMM/MM  Code 
(cpu  secs) 

0° 

353.8 

20/20 

95.7 

5° 

600 

357.3 

50/30 

170.0 

45° 

600 

353.6 

100/20 

334.3 

oo 

CJ\ 

o 

600 

356.8 

200/20 

1980.7 

Table  2.  Wire  Array  in  Slab  Reference  MM  Code  versus  Hybrid  PMM/MM  Code  Time  Comparison 


Figures  52  through  55  show  the  induced  currents  for  a  plane  wave  incident  45°  off  normal. 
Examining  the  figures  reveal  100  and  20  edge  elements  were  sufficient  to  accurately  model  the  leading 
and  trailing  edges  in  this  scattering  problem. 

Finally,  Figures  56  through  59  show  the  currents  induced  by  a  grazing  illumination  of  85°  off 
normal.  In  this  case,  the  author  used  200  and  20  edge  elements  for  the  two  semi-infinite  subproblems. 
Notice  the  code  calculates  the  trailing  edge  elements  accurately.  The  dielectric  slab  decreases  the 
angle  of  incidence  in  the  slab,  thus  avoiding  the  numerical  problems  encountered  in  the  homogeneous 
media  geometry. 

As  for  the  homogeneous  case.  Table  2  contains  the  CPU  times  for  the  reference  MM  and 
hybrid  PMM/MM  based  codes  to  compute  each  scattering  scenario.  Unlike  the  homogeneous  case, 
these  values  represent  a  true  “apples-to-apples”  comparison  since  the  reference  code  used  the  core 
of  the  hybrid  code  with  the  PBF  “turned  off.”  While  these  computation  times  are  not  as  dramatic 
as  the  homogeneous  case,  the  fact  remains  as  the  problem  size  increases,  the  MM  based  system 
sizes  quickly  become  prohibitively  large,  where  as  the  hybrid  PMM/MM  based  system  size  remains 
constant.  Note,  the  large  hybrid  technique  computation  time  for  an  85°  incident  angle  is  a  result  of 
numerical  difficulties  in  computing  the  continuous  spectral  integral  and  not  truly  indicative  of  the 
technique’s  performance.  Asymptotic  evaluation  would  probably  reduce  this  time. 
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Figure  46.  Induced  Current  Magnitude  on  Left  50  Elements  of  Wire  Array  in  Slab  with  Normal 
Illumination 
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Figure  47.  Induced  Current  Phase  on  Left  50  Elements  of  Wire  Array  in  Slab  with  Normal 
Illumination 


83 


Element  Number 


Figure  50.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  in  Slab  with  5°  off 
Normal  Illumination 


Figure  5 1 .  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  in  Slab  with  5°  off  Normal 
Illumination 
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Figure  54.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  in  Slab  with  45°  off 
Normal  Illumination 


Figure  55.  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  in  Slab  with  45°  off  Normal 
Illumination 
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Figure  58.  Induced  Current  Magnitude  on  Right  50  Elements  of  Wire  Array  in  Slab  with  85°  off 
Normal  Illumination 
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Figure  59.  Induced  Current  Phase  on  Right  50  Elements  of  Wire  Array  in  Slab  with  85°  off  Normal 
Illumination 


89 


4.13  Slot  Array  in  a  Homogeneous  Media.  Having  examined  the  two-dimensional  wire 
array,  the  chapter  now  focuses  on  several  slot  array  geometries.  Slotted  ground  planes  are  a  common 
form  of  radome.  Thus,  the  following  plots  demonstrate  the  hybrid  PMM/MM  technique’s  applicability 
to  practical  radome  design  problems. 

As  for  the  wire  arrays,  the  author  validated  the  hybrid  code  against  an  ASM  based  MM  code. 
In  this  case,  two  independent  reference  codes,  one  written  by  the  author  and  one  by  Skinner  [33], 
provide  additional  confidence  in  the  solution  accuracy.  The  reference  code  solutions  represent  the 
“truth”  standard  one  uses  to  judge  the  hybrid  code’s  validity.  One  desires  both  agreement  between 
the  reference  codes  and  agreement  with  the  hybrid  code.  In  terms  of  testing  the  validity  of  the  this 
research,  the  latter  comparison  is  by  far  the  more  important.  Thus,  one  must  examine  the  following 
figures  with  two  separate  questions  in  mind.  One,  do  the  two  reference  codes  agree  reasonably  well? 
Two,  does  the  hybrid  code  accurately  represent  the  reference  code  solution? 

Figure  23,  at  the  beginning  of  Chapter  III,  depicts  the  general  features  of  the  slot  arrays 
examined  in  the  following  sections.  This  particular  section  begins  by  examining  a  100  column  slot 
array  in  a  homogeneous  media.  The  elements  are  z  oriented  slots  0.881  cm  long  and  0.038  cm  wide. 
They  form  slot  columns  with  inner-element  spacing  1.27  cm  in  the  z  direction  and  1.1  cm  in  the  ® 
direction.  The  columns  create  a  non-skewed  grid.  A  15  GHz  plane  wave  illuminates  the  array  at  5°, 
45°,  and  85°  off  normal.  Element  number  one  always  belongs  to  the  leading  edge  of  the  array.  In  all 
cases,  the  author  uses  a  single  PWS  mode  per  slot. 

Figures  60  and  61  show  the  magnitude  and  phase  of  the  currents  induced  on  the  column 
reference  elements.  The  reference  values  are  denoted  by  Xs  and  *s,  and  the  hybrid  PMM/MM  values 
are  denoted  by  Os.  The  array  is  illuminated  by  a  plane  wave  incident  5°  off  normal.  The  hybrid 
technique  used  20  edge  elements  for  both  left  and  right  hand  problems.  First,  note  the  two  reference 
codes  agree  reasonably  well.  The  current  perturbations  near  the  edges  are  virtually  identical  in  the 
two  solutions.  The  main  difference  is  a  slight  offset  between  the  curves.  This  offset  can  be  attributed 
to  the  different  numerical  approaches  the  two  authors  used.  In  fact,  the  good  agreement  is  a  testament 
to  the  validity  of  the  underlying  ASM  approach.  Second,  compare  the  hybrid  solution  to  the  reference 
solutions.  The  agreement  is  excellent.  Since  the  author’s  reference  solution  is  generated  by  the  hybrid 
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code  with  the  PBF  turned  off,  it  provides  a  more  accurate  standard  against  which  to  judge  the  hybrid 
solution. 

Figures  62  and  63  show  the  currents  induced  on  the  same  array  illuminated  by  a  plane  wave 
incident  45°  off  normal.  Here,  the  hybrid  technique  used  25  edge  elements  for  both  left  and  right 
hand  problems.  In  this  case,  the  two  reference  codes  provide  virtually  identical  results.  The  hybrid 
solution  accurately  captures  the  edge  element  behavior  for  those  elements  modeled.  However,  a  few 
more  edge  elements  would  give  a  slightly  more  accurate  PBF  value  for  the  array’s  central  portion. 

Finally,  Figures  64  and  65  show  the  currents  induced  on  the  same  array  illuminated  by  a  plane 
wave  incident  85°  off  normal.  The  hybrid  technique  used  30  edge  elements  for  both  left  and  right  hand 
problems.  The  hybrid  solution  is  a  good  approximation  to  the  reference  codes’  solutions.  However, 
as  for  the  wire  array,  the  author’s  code  has  difficulties  capturing  the  trailing  edge  behavior  at  extreme 
grazing  angles.  Fortunately,  in  the  homogeneous  media  where  the  problem  occurs,  these  currents 
always  have  a  much  smaller  magnitude  than  the  leading  edge  currents,  mitigating  their  influence  in 
the  far-held.  Figure  66  demonstrates  this  fact  with  a  monostatic  radar  echo  area  (REA)  plot  over  the 
angle  range  from  0°  to  88°.  Although  not  implemented  in  the  author’s  codes,  it  is  apparent  from  the 
hgure  one  could  save  computation  time  with  little  sacrihce  in  accuracy  by  only  computing  the  leading 
edge  semi-infinite  problem  at  near  grazing  incident  angles. 
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Figure  60.  Induced  Current  Magnitude  on  100  Element  Slot  Array  with  5°  off  Normal  Illumination 


Figure  61.  Induced  Current  Phase  on  100  Element  Slot  Array  with  5°  off  Normal  Illumination 
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Figure  62.  Induced  Current  Magnitude  on  1 00  Element  Slot  Array  with  45°  off  Normal  Illumination 


Figure  63.  Induced  Current  Phase  on  100  Element  Slot  Array  with  45°  off  Normal  Illumination 
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Figure  64.  Induced  Current  Magnitude  on  100  Element  Slot  Array  with  85°  off  Normal  Illumination 
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Figure  65.  Induced  Current  Phase  on  100  Element  Slot  Array  with  85°  off  Normal  Illumination 


94 


Monostatic  Angle  Eta  (Degrees) 


Figure  66.  REA  pattern  at  15  GHz  for  100  Column  Slot  Array  in  Homogeneous  Media 
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4.1.4  Slot  Array  in  a  Stratified  Media.  This  section  adds  a  symmetric  dielectric  sandwich 
around  the  the  array  examined  in  the  previous  section.  The  top  and  bottom  layers  are  0.55875  cm 
thick  and  have  relative  permittivities  of  2.0.  Once  again,  the  differences  between  the  reference  codes 
can  be  attributed  to  the  authors’  numerical  implementations.  Specifically,  the  authors  use  different 
approaches  to  handling  the  surface  wave  singularities  introduced  by  the  dielectric  layers.  For  the 
single  layer  case,  Skinner  uses  an  analytic  extraction  technique  to  remove  the  T-factor  singularities. 
This  contrasts  with  the  lossy  dielectric  approach  described  earlier.  It  is  difficult  to  ascertain  which 
solution  is  closer  to  the  truth,  but  again,  the  main  issue  is  how  well  the  hybrid  technique  captures  the 
MM  solution.  On  this  matter,  the  following  figures  speak  for  themselves. 

Figures  67  and  68  show  the  current  magnitude  and  phase  induced  on  the  column  reference 
elements  for  a  5°  off  normal  incident  angle.  The  hybrid  technique  used  20  edge  elements  for  both  left 
and  right  hand  problems.  As  for  the  homogeneous  case,  the  hybrid  PMM/MM  accurately  captures 
the  induced  currents,  as  differences  in  Figure  67  have  been  magnified  by  an  expanded  scale. 

Figures  69  and  70  show  the  current  magnitude  and  phase  induced  on  the  column  reference 
elements  for  a  45°  off  normal  incident  angle.  The  hybrid  technique  used  20  edge  elements  for  both 
left  and  right  hand  problems.  Notice  how  the  dielectric  sandwich  has  modified  the  trailing  edge 
current  behavior.  Once  again,  the  hybrid  technique  compares  well  with  the  reference  solution. 

Figures  71  and  72  show  the  current  magnitude  and  phase  induced  on  the  column  reference 
elements  for  a  85°  off  normal  incident  angle.  The  hybrid  technique  used  20  edge  elements  for  both 
left  and  right  hand  problems.  Clearly,  the  leading  edge  requires  more  edge  elements.  Note  how 
the  PBF  magnitudes  for  the  two  half  problems  differ.  This  provides  a  good  indication  more  edge 
elements  are  necessary  to  improve  the  solution  accuracy.  In  addition,  one  may  note  the  trailing  edge 
difficulties  associated  with  the  homogeneous  array  are  virtually  eliminated  by  the  dielectric  cladding, 
as  seen  earlier  with  the  wire  array.  This  can  be  attributed  to  refraction  as  the  plane  wave  is  bent 
toward  normal  upon  entering  the  dielectric  slab. 
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Figure  67.  Induced  Current  Magnitude  on  100  Element  Slot  Array  in  Stratified  Media  with  5°  off 
Normal  Illumination 


Figure  68.  Induced  Current  Phase  on  100  Element  Slot  Array  in  Stratified  Media  with  5°  off  Normal 
Illumination 
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Figure  69.  Induced  Current  Magnitude  on  100  Element  Slot  Array  in  Stratified  Media  with  45°  off 
Normal  Illumination 


Figure  70.  Induced  Current  Phase  on  100  Element  Slot  Array  in  Stratified  Media  with  45°  off 
Normal  Illumination 


Figure  71. 


Figure  72. 


Induced  Current  Magnitude  on  100  Element  Slot  Array  in  Stratified  Media  with  85°  off 
Normal  Illumination 


Induced  Current  Phase  on  100  Element  Slot  Array  in  Stratified  Media  with  85°  off 
Normal  Illumination 
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4.1.5  Slot  Array  with  Loaded  Edge  Elements.  The  ability  to  modify  the  edge  element 
currents  via  loading  is  an  important  tool  to  the  radome  designer.  This  section  presents  the  results  of 
applying  r-card  to  the  first  and  last  nine  columns  of  the  homogeneous  media  slot  array  described  in 
section  4.1.3.  Specifically,  the  outer  three  edge  element  modes  are  loaded  with  231  g  r-card,  the  next 
three  have  6955g  r-card,  and  the  final  three  have  139100  r-card. 

Figures  73  through  78  show  the  effect  of  this  treatment  on  the  array,  once  again  illuminated  at 
5°,  45°,  and  85°  off  normal.  Note  the  two  reference  solutions  and  the  hybrid  solution  all  agree  fairly 
well,  except  for  the  aforementioned  difficulty  with  the  trailing  edge.  The  slight  discrepancies  in  the 
loaded  element  currents  are  a  result  of  Skinner  approximating  the  slot  current  mode  as  a  constant  to 
calculate  the  lumped  load.  This  contrasts  with  the  author’s  method  described  in  section  4.1.5.  Since 
Skiimer  uses  a  different  approach  to  determine  the  proper  lumped  load  to  apply  to  the  admittance 
matrix,  it  is  encouraging  to  see  the  same  end  result.  This  example  demonstrates  the  code’s  usefulness 
in  designing  array  edge  treatments. 

Section  4.1  validated  the  hybrid  PMM/MM  technique  by  comparing  the  author’s  numerical 
implementation  to  other  reference  codes  based  on  the  MM.  Since  the  reference  data  is  calculated 
data,  the  actual  induced  array  currents  can  be  compared,  providing  the  most  rigorous  test  of  the 
hybrid  technique’s  accuracy.  The  results  clearly  demonstrate  the  technique’s  validity  and  usefulness 
in  predicting  large  finite  array  performance.  To  add  another  layer  of  confidence,  the  author  also 
measured  several  radomes  and  compared  the  measured  echo  area  with  that  calculated  by  the  hybrid 
PMM/MM  code.  These  results  are  presented  in  the  next  section. 
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Figure  73. 


Figure  74. 


Induced  Current  Magnitude  on  100  Element  Loaded  Slot  Array  with  5°  off  Normal 
Illumination 


Induced  Current  Phase  on  100  Element  Loaded  Slot  Array  with  5°  off  Normal 
Illumination 


101 


8 


7h 


■ 

■ 

■ 

■ 


^4- 

1 

3- 


X  Reference  MM  Solution  (Skinner) 
*  Reference  MM  Solution  (Collins) 
o  Hybrid  PMM/MM  Solution 


lh 


m 


10  20  30  40  50  60  70  80  90  100 

Element  Number 


Figure  75.  Induced  Current  Magnitude  on  100  Element  Loaded  Slot  Array  with  45°  off  Normal 
Illumination 


Figure  76.  Induced  Current  Phase  on  100  Element  Loaded  Slot  Array  with  45°  off  Normal 
Illumination 
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Figure  77. 


Figure  78. 
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Induced  Current  Magnitude  on  100  Element  Loaded  Slot  Array  with  85°  off  Normal 
Illumination 


Induced  Current  Phase  on  100  Element  Loaded  Slot  Array  with  85°  off  Normal 
Illumination 
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4.2  Measurement  Comparison 

Recall  the  geometries  of  concern  consist  of  finite  planar  arrays  made  from  slot  colunms 
infinite  in  one  dimension.  In  addition,  the  slot  arrays  exist  on  infinite  ground  planes.  These  facts 
pose  a  somewhat  difficult  measurement  problem,  since  infinite  structures  are  hard  to  fit  in  a  finite 
measurement  facility.  The  challenge,  then,  is  to  approximate  these  structures  with  the  finite  resources 
available.  The  author  attempted  to  design  finite  radomes  that  were  large  enough  in  the  i  direction  to 
approximate  infinite  columns.  Fortunately,  the  thin  slot  pattern  factors  dramatically  reduce  the  edge 
effects  in  the  £  direction,  allowing  a  relatively  short  column  to  adequately  approximate  an  infinite 
column. 

The  Mission  Research  Corporation  of  Dayton,  Ohio  manufactured  the  two  different  radomes 
the  author  used  during  the  measurement  validation  phase.  The  radomes  were  measured  at  the  Wright 
Laboratory  Multi-Spectral  Measurement  Facility,  Wright-Patterson  AFB,  Ohio.  Wright  Laboratory 
also  provided  the  test  fixture  used  during  the  measurements. 

The  radomes  consist  of  skewed  grids  of  z  oriented  slots  0.881  cm  long  and  0.038  cm  wide. 
The  slots  form  columns  with  inner-element  spacing  1.27  cm  in  the  z  direction  and  1.1  cm  in  the  ® 
direction.  Each  radome  has  42  slot  columns,  creating  a  window  approximately  45  cm  wide  (x)  by 
51  cm  tall  (z).  The  grids  have  a  45°  skew  angle.  The  first  radome  is  simply  a  thin  metal  screen.  The 
second  radome  is  backed  by  a  0.15875  cm  thick  glass-epoxy  layer  with  a  relative  permittivity  of  4.5. 
The  author  measured  each  radome  with  and  without  a  similar  r-card  edge  treatment  to  that  described 
in  the  previous  section.  In  this  case,  resistivities  of  500g,  750§,  and  1500g  were  used  for  a  total  of 
four  radome  configurations. 

Figure  79  shows  Wright  Laboratory’s  test  fixture.  The  test  fixture  attempts  to  simulate  an 
infinite  ground  plane  by  mitigating  any  backscatter  due  to  various  mechanisms  such  as  tip  scatter  and 
creeping  waves.  Its  size  also  aids  in  gating  out  any  interactions  between  the  radome  and  test  body. 

Figures  82  and  83  show  a  radome  mounted  in  the  test  body  and  the  cavity  beneath  the  radome, 
respectively.  Since  the  radomes  are  transmissive  within  the  band  of  interest  (see  Figures  80  and 
81),  the  cavity  creates  a  contamination  source  by  re-radiating  energy  coupled  into  the  cavity  back 
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Figure  79.  Radome  Measurement  Test  Fixture 


toward  the  radar.  To  mitigate  this  problem,  the  author  filled  the  cavity  with  RAM  to  absorb  most  of 
the  radiation  transmitted  into  the  cavity.  As  can  be  seen  in  the  following  plots,  this  technique  was 
generally  quite  effective. 

The  final  photograph  shows  a  typical  r-card  treatment.  The  r-card  consists  of  a  thin  polymer 
loaded  with  resistive  material  such  as  carbon.  As  shown  in  Figure  84,  one  attaches  the  film  to  the 
radome  with  spray  adhesive.  One  then  uses  a  putty  knife  to  smooth  out  any  bubbles  trapped  by  the 
film.  The  author  found  the  r-card  material  easy  to  work  with  and  very  consistent  from  measurement 
to  measurement. 
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Figure  81.  Dielectric  Backed  Radome  Transmission  for  Three  Angles  of  Incidence 
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Figure  82.  Radome  Mounted  in  Test  Fixture 


Figure  83.  Test  Fixture  Cavity  Treatment 
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Radome  Configuration 

Frequency  Scans  (GHz) 

Patterns 

Metal 

Metal  with  r-card 
Backed  Metal 
Backed  Metal  with  r-card 

2  to  18(0.1) 

2  to  18  (0.1) 

2  to  18  (0.1) 

2  to  18  (0.1) 

-90°  to  -30°  (0.2°) 
-90°  to  -30°  (0.2°) 
-90°  to  -30°  (0.2°) 
-90°  to  -30°  (0.2°) 

Table  3.  Radome  Measurement  Matrix 

The  Multi-Spectral  Measurement  Facility  is  able  to  take  a  tremendous  amount  of  data  iu  a  very 
short  time.  Table  3  contains  the  test  matrix  for  the  four  radome  configurations.  Frequency  versus 
angle  plots  provide  a  convenient  way  to  view  a  large  amount  of  data.  Figures  85  through  88  show 
the  radar  echo  area  (REA)  for  the  four  radome  configurations.  Note  that  0°  elevation  corresponds  to 
a  grazing  angle  of  incidence.  The  author  chose  to  focus  on  the  three  frequencies,  13  GHz,  15  GHz, 
and  17  GHz,  to  validate  the  hybrid  PMM/MM  technique.  The  remainder  of  this  section  examines 
each  of  the  four  radomes  at  these  frequencies. 
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Figure  85.  Frequency  versus  Angle  Plot  for  Radome  Configuration  Number  1 


Figure  86.  Frequency  versus  Angle  Plot  for  Radome  Configuration  Number  2 
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Figure  87.  Frequency  versus  Angle  Plot  for  Radome  Configuration  Number  3 


Figure  88.  Frequency  versus  Angle  Plot  for  Radome  Configuration  Number  4 
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Figure  89.  REA  pattern  at  13  GHz  for  Metal  Radome  in  Stratified  Media 

4.2.1  Dielectric  Backed  Metal  Radome.  Figures  89  through  91  depict  monostatic  REA 
patterns  at  13  GHz,  15  GHz,  and  17  GHz  respectively.  Each  plot  contains  measured  data,  reference 
code  data  calculated  using  the  MM,  and  hybrid  PMM/MM  data  based  on  10  edge  columns  on  either 
side  of  the  array.  Several  features  are  worth  noting  in  this  figure  series.  First  of  all,  the  two  reference 
codes  agree  quite  well  in  all  the  figures.  While  not  perfect,  they  also  do  a  good  job  capturing  the 
behavior  of  measured  data,  improving  with  increasing  frequency.  Second,  the  hybrid  PMM/MM  data 
agrees  with  the  reference  data  fairly  well  over  some  regions  but  misses  the  mark  in  others.  As  an 
example,  notice  how  the  hybrid  code  misses  the  lobe  at  —  50°  in  Figure  90  completely.  To  understand 
what  is  happening  here,  one  must  examine  the  induced  scattering  currents.  The  following  current 
plots  highlight  one  of  the  code’s  useful  features.  Namely,  the  ability  to  access  the  scattering  currents. 
This  ability  gives  the  radome  designer  insight  into  what  is  producing  a  particular  REA  pattern  or 
perhaps  the  effectiveness  of  a  particular  edge  treatment. 

Figures  92  and  93  present  the  induced  current  magnitudes  calculated  by  the  author’s  MM  and 
hybrid  PMM/MM  codes  at  15  GHz.  In  the  MM  solution,  notice  how  the  dielectric  layer  creates  a 
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Figure  90.  REA  pattern  at  15  GHz  for  Metal  Radome  in  Stratified  Media 

strong  coupling  between  the  slot  columns.  This  is  evidenced  by  the  absence  of  a  central  Floquet  region, 
even  at  angles  far  away  from  grazing.  The  hybrid  technique  would  be  hard  pressed  to  accurately 
model  the  central  currents  for  a  radome  this  narrow.  Even  more  noticeable  is  the  missing  ridge  at 
around  —50°  in  the  hybrid  solution.  The  author  traced  the  problem  to  the  continuous  spectral  integral 
evaluations  for  the  PBF.  In  this  case,  the  dielectric  layer  produces  T-factor  poles  which  lie  either  near 
grating  lobe  poles  or  near  the  branch  points  in  the  complex  integration  plane.  The  assumptions  the 
author  used  to  analytically  extract  the  grating  lobe  pole  may  be  violated,  necessitating  an  alternative 
approach  to  dealing  with  the  singularity. 
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Figure  91. 


Figure  92.  Induced  Currents  on  Metal  Radome  in  Stratified  Media  at  15  GHz 


Figure  93.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  Metal  Radome  in  Stratified 
Media  at  15  GHz 
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4.2.2  Dielectric  Backed  Metal  Radome  with  R-Card  Edge  Treatment.  Figures  94  through 
96  depict  the  r-card  treated  dielectric  backed  metal  radome ’s  monostatic  REA  patterns  at  13  GHz, 
1 5  GHz,  and  1 7  GHz  respectively.  As  mentioned  earlier,  the  hybrid  PMM/MM  code  ran  into  numerical 
difficulties.  Unfortunately,  Skinner’s  reference  MM  code  also  had  problems.  Figure  94  contains  the 
only  usable  data,  which  may  be  suspect.  At  the  higher  frequencies,  Skinner’s  code  returned  NaNs 
(Not  a  Number),  indicating  numerical  instability.  Despite  the  lack  of  computed  reference  data,  the 
author’s  MM  code  agrees  quite  well  with  the  measured  data.  One  may  note  how  the  edge  treatments 
have  widened  the  lobing  structure,  while  reducing  the  side  lobes.  This  behavior  can  be  accounted  for 
by  examining  the  currents  in  Figures  97  and  98.  The  r-card  has  partially  shorted  out  the  outer  edge 
currents,  effectively  reducing  the  width  of  the  radome. 


Figure  94.  REA  pattern  at  13  GHz  for  R-Card  Treated  Metal  Radome  in  Stratified  Media 


116 


Monostatic  Angle  Eta  (Degrees) 


Figure  96.  REA  pattern  at  17  GHz  for  R-Card  Treated  Metal  Radome  in  Stratified  Media 
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Figure  97.  Induced  Currents  on  R-Card  Treated  Metal  Radome  in  Stratified  Media  at  15  GHz 
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Figure  98.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  R-Card  Treated  Metal  Radome 
in  Stratified  Media  at  15  GHz 


4.2.3  Metal  Radome  in  Homogeneous  Media.  Figures  99  through  101  contain  monostatic 
REA  patterns  for  the  metal  radome.  The  frequencies  of  interest  are  again  13  GHz,  15  GHz,  and 
17  GHz.  Each  plot  contains  measured  data,  reference  code  data,  and  hybrid  PMM/MM  data  based 
on  16  edge  columns  on  either  side  of  the  array.  In  examining  the  measured  data,  there  appears  to 
be  a  greater  discrepancy  with  the  predicted  results  than  observed  with  the  dielectric  backed  radome. 
Unfortunately,  the  measured  data  for  the  metal  radome  is  somewhat  suspect.  There  are  several 
possible  explanations.  First,  the  metal  foil’s  flimsiness  may  have  contributed.  In  order  to  take 
background  measurements  for  background  subtraction,  the  author  lightly  adhered  aluminum  foil 
to  the  metal  radome  surface.  Removing  the  foil  without  disturbing  the  metal  radome  proved  nearly 
impossible.  As  a  result,  the  background  subtraction  introduced  phase  errors  which  may  have  corrupted 
the  measurements.  Second,  the  test  fixture  may  have  been  tilted  slightly  away  from  perpendicular  to 
the  radar’s  line-of-sight.  The  effect  would  be  a  reduction  in  peak  height.  Third  and  most  probable, 
the  cavity  return  may  have  corrupted  the  measurement.  Since  the  metal  radome  is  more  transmissive 
in  the  frequency  band  of  interest,  the  RAM  inside  the  cavity  wasn’t  able  to  completely  absorb  the 
coupled  radiation  causing,  the  gross  lobing  structure  in  Figure  99  and  the  grating  lobe  peak  height 
discrepancies  in  Figures  100  and  101. 

Turning  now  to  the  predicted  results,  note  once  again,  the  reference  codes  produce  virtually 
identical  solutions.  One  may  also  note  the  excellent  agreement  with  the  hybrid  PMM/MM  predictions. 
The  only  areas  of  deviation  are  near  grazing.  Based  on  Figures  102  through  107,  it  is  clear  the  problem 
at  grazing  is  caused  once  again  by  numerical  problems  at  the  trailing  edge.  Also,  as  the  incident 
angle  sweeps  off  normal,  the  Floquet  portion  at  the  center  of  the  array  shrinks  to  the  point  where  the 
PBF  is  no  longer  useful.  This  effect  is  evidenced  in  the  REA  pattern  by  the  hybrid  solution  departing 
from  the  MM  solution  as  the  incident  angle  approaches  grazing.  Before  leaving  these  figures,  it  is 
interesting  to  note  the  current’s  behavior  as  the  grating  lobe  enters  the  visible  region  near  13  GHz  and 
moves  toward  normal  with  increasing  frequency.  The  currents  on  one  side  of  the  grading  lobe  have 
a  dramatically  different  behavior  than  those  on  the  other  side. 
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Figure  99.  REA  pattern  at  13  GHz  for  Metal  Radome  in  Homogeneous  Media 
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Figure  100.  REA  pattern  at  15  GHz  for  Metal  Radome  in  Homogeneous  Media 
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Figure  101.  REA  pattern  at  17  GHz  for  Metal  Radome  in  Homogeneous  Media 
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Figure  102.  Induced  Currents  on  Metal  Radome  in  Homogeneous  Media  at  13  GHz 
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Figure  103.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  Metal  Radome  in  Homoge 
neous  Media  at  13  GHz 
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Figure  104.  Induced  Currents  on  Metal  Radome  in  Homogeneous  Media  at  15  GHz 
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Figure  105.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  Metal  Radome  in  Homoge 
neous  Media  at  15  GHz 
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Figure  107.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  Metal  Radome  in  Homoge¬ 
neous  Media  at  17  GHz 
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4.2.4  Metal  Radome  in  Homogeneous  Media  with  R-Card  Edge  Treatment.  Figures  108 
through  110  depict  the  r-card  treated  metal  radome ’s  monostatic  REA  patterns  at  13  GHz,  15  GHz, 
and  17  GHz  respectively.  Each  plot  contains  measured  data,  reference  code  data  calculated  using  the 
MM,  and  hybrid  PMM/MM  data  based  on  16  edge  columns  on  either  side  of  the  array.  Once  again, 
the  calculated  data  does  a  good  job  predicting  the  radome’s  general  performance.  The  measured  data 
may  have  errors  due  to  measurement  difficulties  as  evidenced  by  the  grating  lobe  peak  discrepancies. 
Finally,  these  figures  show  where  the  hybrid  PMM/MM  solution  best  represents  the  actual  array 
currents.  Not  surprisingly,  this  occurs  when  the  currents  have  a  good  central  Floquet  region. 


Figure  108.  REA  pattern  at  13  GHz  for  R-Card  Treated  Metal  Radome  in  Homogeneous  Media 
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Figure  109.  REA  pattern  at  15  GHz  for  R-Card  Treated  Metal  Radome  in  Homogeneous  Media 
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Figure  110.  REA  pattern  at  17  GHz  for  R-Card  Treated  Metal  Radome  in  Homogeneous  Media 
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Figure  111.  Induced  Currents  on  R-Card  Treated  Metal  Radome  in  Homogeneous  Media  at  13  GHz 
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Figure  1 12.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  R-Card  Treated  Metal 
Radome  in  Homogeneous  Media  at  13  GHz 
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Figure  113.  Induced  Currents  on  R-Card  Treated  Metal  Radome  in  Homogeneous  Media  at  15  GHz 
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Figure  114.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  R-Card  Treated  Metal 
Radome  in  Homogeneous  Media  at  15  GHz 
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Figure  1 15.  Induced  Currents  on  R-Card  Treated  Metal  Radome  in  Homogeneous  Media  at  17  GHz 


Figure  116.  Hybrid  PMM/MM  Approximation  to  Induced  Currents  on  R-Card  Treated  Metal 
Radome  in  Homogeneous  Media  at  17  GHz 
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43  Array  Size  Scaling 

In  this  final  section,  the  author  presents  an  example  array  problem  demonstrating  the  effect 
of  array  size  on  computation  time  for  both  MM  and  hybrid  MM/PMM  techniques.  It  addresses  the 
question  of  why  one  would  be  interested  using  the  hybrid  technique  if  the  MM  analysis  provides 
adequate  solutions.  To  provide  a  fair  comparison,  the  author  used  the  same  code  with  the  PBF  turned 
on  and  off.  This  approach  eliminated  any  differences  due  to  the  MM’s  numerical  implementation. 
The  author  used  the  same  slotted  array  parameters  as  those  in  the  homogeneous  media  examples 
presented  earlier.  The  array  is  again  illuminated  by  a  15  GHz  plane  wave,  5°  off  normal.  The  only 
variable  is  the  number  of  slot  columns  comprising  each  array. 

Table  4  presents  the  comparison  between  the  two  techniques  for  radomes  with  100  to  900 
slot  columns.  The  codes  were  run  on  a  Sun  Microsystems  Sparc20  Workstation.  In  obtaining  the 
data  for  the  table,  the  author  assumed  twenty  edge  elements  were  sufficient  to  capture  all  the  current 
perturbations  due  to  the  array  edges.  Figure  60  clearly  demonstrates  this  assumption  is  valid.  Notice 
that  while  the  MM  is  faster  than  the  hybrid  PMM/MM  for  a  100  column  array,  the  hybrid  technique 
quickly  surpasses  the  MM  for  larger  arrays.  One  may  also  observe  that,  while  the  MM  matrix  fill  time 
claims  the  lion’s  share  of  the  computation  time  early  on,  the  system  solve  time  quickly  becomes  the 
dominant  bottleneck.  Based  on  the  table’s  data,  the  solve  time  increases  as  approximately  where 
N  indicates  the  system  size.  This  is  in  contrast  to  the  hybrid  PMM/MM  technique,  whose  computation 
time  scales  as  approximately  N^-^.  This  study  clearly  demonstrates  one  of  the  technique’s  primary 
advantages  over  conventional  approaches.  Namely,  the  hybrid  PMM/MM  technique  is  able  to  side¬ 
step  the  link  between  array  size  and  system  size  by  lumping  the  array’s  Floquet  elements  into  a 
single  unknown  element.  As  stated  earlier  for  the  wire  array,  the  computational  improvement  for 
calculations  involving  many  monostatic  angles  may  not  be  as  dramatic  since  the  hybrid  method 
requires  recalculation  of  one  row  in  the  coupling  matrix  for  every  incident  angle. 

This  chapter  represents  the  culmination  of  the  previous  chapters  by  putting  the  theory  into 
practice.  The  previous  data  sets  clearly  demonstrate  where  the  hybrid  PMM/MM  technique  is  useful 
as  well  as  where  it’s  use  is  inappropriate.  The  overall  agreement  between  the  reference  code  and 
measured  data  give  a  high  confidence  level  in  the  technique’s  validity.  In  general,  where  ever  an 
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MM  Code 

MM  Code 

MM  Code 

MM  Code 

Unknowns 

Fill  (cpu  secs) 

Solve  (cpu  secs) 

Total  (cpu  secs) 

6.00 

0.18 

6.27 

12.07 

1.32 

13.48 

300 

18.65 

4.30 

23.05 

400 

24.92 

9.98 

35.02 

500 

30.50 

19.73 

50.38 

600 

36.70 

34.05 

70.95 

700 

54.58 

54.28 

109.12 

800 

49.58 

81.35 

131.12 

900 

56.97 

116.02 

173.25 

PMM/MM  Code 

PMM/MM  Code 

PMM/MM  Code 

PMM/MM  Code 

Unknowns:  Left/Right 

Fill  (cpu  secs) 

Solve  (cpu  secs) 

Total  (cpu  secs) 

4.82/4.78 

0.00/0.02 

9.67 

4.80/4.78 

0.02/0.00 

9.68 

20/20 

4.87/4.85 

0.00/0.02 

9.85 

20/20 

4.80/4.80 

0.00/0.00 

9.72 

20/20 

4.80/4.77 

0.02/0.00 

9.73 

20/20 

4.82/4.80 

0.00/0.02 

9.77 

20/20 

4.85/4.85 

0.02/0.00 

9.92 

20/20 

4.82/4.87 

0.03/0.02 

9.93 

20/20 

4.78/4.80 

0.00/0.02 

9.82 

Table  4.  Slot  Array  in  Homogeneous  Media  Reference  MM  Code  versus  Hybrid  PMM/MM  Code 
Time  Comparison 


array  is  wide  enough  to  have  a  well  developed  Floquet  region,  the  hybrid  PMM/MM  technique  is 
computationally  superior  to  the  MM  formulation. 
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V.  Conclusion 


Of  making  many  books  there  is  no  end,  and  much  smdy  wearies  the  body.  Now  all  has 
been  heard;  here  is  the  conclusion  of  the  matter:  Fear  God  and  keep  His  conunandments, 
for  this  is  the  whole  duty  of  man. 

-  King  Solomon 

This  dissertation  develops  a  hybrid  periodic  moment  method/moment  method  (PMM/MM) 
based  solution  for  the  electromagnetic  scattering  from  semi-infinite  arrays  of  wires  and  slots,  possibly 
embedded  in  a  stratified  dielectric  media.  These  geometries  are  of  interest  in  developing  treatments 
to  control  the  edge  diffraction  from  large  finite  arrays.  While  a  growing  body  of  finite  array  research 
exists,  none  incorporates  the  physical  basis  function  (PBF)  concept  presented  here.  The  dissertation 
clearly  demonstrates  the  benefits  of  the  hybrid  approach.  With  the  exception  of  a  few  instances  where 
the  author’s  numerical  implementation,  not  the  formulation  itself,  produces  inaccurate  results,  the 
hybrid  code  accurately  generates  the  array’s  scattering  currents  in  much  less  time  than  the  reference 
moment  method  (MM)  codes. 

In  the  previous  pages,  the  author  develops  each  component  of  a  complete  moment  method  (MM) 
solution  for  both  wire  and  slot  arrays.  Examining  each  piece  sparately,  the  author  fills  the  coupling 
matrix  consisting  of  individual  edge  element  basis  functions  and  a  PBF  capturing  the  Floquet  behavior 
away  from  the  array’s  edge.  With  an  eye  toward  an  efficient  numerical  implementation,  both  spatial 
and  spectral  domain  forms  are  presented,  along  with  various  summation  acceleration  techniques. 
Dielectric  boundaries  are  accounted  for  via  Fresnel  reflection  and  transmission  coefficients.  The 
author  then  turns  to  the  excitation  vector  entries,  again  accounting  for  the  dielectric  boundaries  in 
the  spectral  domain  via  Fresnel  coefficients.  Solving  the  system,  the  author  discusses  the  scattering 
currents  and  their  relationship  to  the  scattered  far-fields  and  echo  widths.  Various  related  ideas  are  also 
discussed  including  PBF  testing  and  use  in  the  excitation  vector,  a  priori  determination  of  the  required 
number  of  edge  elements,  and  edge  element  loading.  Finally,  the  hybrid  PMM/MM  formulation  is 
validated,  using  both  independent  reference  solutions  and  actual  array  measurements.  The  author 
tests  an  r-card  based  edge  treatment  using  the  hybrid  formulation.  The  resulting  scattering  currents 
and  field  patterns  demonstrate  the  formulation’s  usefulness  in  array  design. 
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This  work  represents  the  first  attempt  to  incorporate  a  periodic  basis  function  with  individual 
edge  basis  functions  in  an  effort  to  reduce  the  computational  burden  associated  with  large  arrays.  The 
computational  savings  demonstrated  here  suggest  several  areas  for  future  research,  taking  advantage 
of  this  promising  technique.  First,  the  hybrid  technique  could  be  applied  to  more  complex  arrays.  The 
author  necessarily  restricted  the  array  complexity  in  order  to  focus  on  the  PBF  implementation  issues. 
The  technique’s  success  motivates  exploring  more  general  element  shapes  and  orientations  in  multi¬ 
array  structures.  The  PBF  could  also  be  used  in  orthogonal  directions,  extending  the  finite-by-infinite 
arrays  examined  here  to  finite-by-finite  geometries.  Second,  the  numerical  difficulties  the  author 
encountered  point  to  explorations  of  alternative  numerical  PBF  implementations.  In  particular,  one 
may  want  to  examine  numerical  integration  algorithms  tailored  to  the  continuous  spectral  integrand’s 
form.  The  asymptotic  analysis  for  the  homogeneous  media  case  presented  in  Appendix  C  could 
be  developed  for  stratified  media  to  further  increase  computational  efficiency.  Additionally,  there 
may  be  different  approaches  to  handling  the  surface  wave  and  grating  lobe  singularities.  Finally,  the 
grazing  angle  problem  caused  by  the  epsilon  summation  acceleration  algorithm  could  be  avoided  by 
implementing  an  alternative  algorithm  such  as  those  discussed  in  Appendix  B.  The  last  research  area 
involves  applying  the  new  hybrid  technique  to  the  edge  treatment  design  problem.  One  could  easily 
take  advantage  of  the  technique’s  speed  by  incorporating  it  in  an  automated  optimization  scheme. 
Regardless  of  the  optimization  approach  (gradient  search,  genetic  algorithm,  etc.),  the  parameter 
space  can  be  searched  much  quicker  due  to  the  faster  minimization  function  calculations.  Of  course, 
a  priori  knowledge  of  the  required  number  of  edge  elements  would  also  help  speed  the  process.  One 
could  incorporate  an  approximation  such  as  the  Gibbsian  model  discussed  at  the  end  of  Chapter  II. 
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Appendix  A.  One-Sided  Poisson  Sum  Formula 


The  coupling  matrix  entries  associated  with  periodic  structures  often  involve  slowly  convergent 
infinite  summations.  The  Poisson  sum  formula  provides  a  convenient  means  to  transform  these 
summations  into  a  more  rapidly  convergent  form.  In  addition,  the  resulting  plane  wave  representation 
allows  the  introduction  of  stratified  dielectric  layers  with  relative  ease. 

The  semi-infinite  array  geometries  of  interest  in  this  research  lead  to  one-sided  summations. 
Unfortunately,  the  original  form  of  the  Poisson  sum  is  inappropriate  in  these  cases.  To  alleviate  this 
difficulty,  a  one-sided  version  of  the  Poisson  sum  formula  is  derived  which  can  be  used  to  transform 
these  one-sided  sums  into  equivalent,  more  rapidly  convergent  plane  wave  representations. 

One  begins  the  derivation  by  considering  the  normal  two-sided  Poisson  sum  formula  given  by 
Papoulis  [29] 

00  00 

Y,  =tY  /(<  ±  kT)  (174) 

3=— 00  k=—oo 

where  f{t)  and  F{u))  denote  a  function  and  its  Fourier  transform  defined  by 
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and  T  =  Note  the  left  hand  side  of  Equation  (174)  contains  a  function  with  the  linear  phase 
progression  characteristic  of  a  spatial  summation,  making  it  an  appropriate  choice  for  conversion  to 
a  one-sided  summation.  This  is  accomplished  by  defining  a  continuous  function  G{uj)  such  that 


F{u;)  =  G(u;)  ■  (/(u;  + 


(177) 
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where  U{uj)  is  the  Heaviside  unit  step  function 


U{u) 


1,  w  >  0 
0,  w  <  0 


(178) 


The  left  hand  side  of  Equation  (174)  can  now  be  written  as 


q=0 


(179) 


which  can  be  identified  with  the  spatial  summations  encountered  in  this  dissertation. 

Now  lets  examine  the  effect  of  this  substitution  on  the  right  hand  side  of  the  double-sided  Poisson 
sum  formula.  One  needs  to  determine  the  inverse  Fourier  transform  of  the  function  G{lj)-U{u}  +  ^). 
Designating  u(t)  and  g{t)  as  the  inverse  Fourier  transforms  of  U{u))  and  G{uj)  respectively,  one  has 
from  the  convolution  property  of  the  Fourier  transform  [35] 


OO 

/ 


G{(jo)  •  U (u)  /  g{T)u{t  —  T)dT 


(180) 


In  addition,  the  frequency  shifting  property  of  the  Fourier  transform  gives  [35] 


U{lj +  '^)  u{t)e  (181) 

Substituting  Equation  (181)  into  Equation  (180) 

OO 

G{bS)  •  U{(jJ  +  •M-  J  g{T)u{t  —  =  f{t)  (182) 

—  OO 


At  this  point,  one  explicitly  introduces  the  inverse  Fourier  transform  of  the  step  function  (Recall 
the  step  function  is  in  the  frequency  domain).  The  duality  property  of  Fourier  transform  pairs  applied 
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to  the  generalized  Fourier  transform  of  the  step  function  results  in  [35] 

/  X  S(t)  1 

u(t)  =  — - — — 

^  ^  2  j27rt 

One  can  substitute  Equation  (183)  in  Equation  (182)  to  obtain 


(183) 


^  2  2%  J  t-T 


based  on  the  sifting  property  of  the  Dirac  delta  function.  The  remaining  integral  must  be  interpreted 
as  a  principal  value  integral  to  maintain  the  one-to-one  correspondence  between  the  Fourier  transform 
pair. 

With  Equation  (184),  one  can  express  the  right  hand  side  in  terms  of  the  new  function  g{t). 
Substituting  Equation  (184)  into  the  original  Poisson  sum  formula  leads  to  the  new  one-sided  Poisson 
sum  formula 


?=0  fc=-oo  {  -oo  ) 


(185) 


which  is  used  throughout  this  research.  Note  the  choice  of  plus  or  minus  signs  is  completely  arbitrary 
in  implementing  this  formula  as  long  as  the  choice  is  applied  consistently. 

The  following  example  illustrates  Equation  (185)’s  implementation.  Consider  the  function, 

,  ,  sin(l) 

9{t)  =  (186) 

and  its  Fourier  transform  [35] 


G{lo)  =  rect((u)  = 


1,  M<1/2 
0,  Iwl  >  1/2 
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Inserting  these  functions  into  Equation  (185),  one  obtains 


oo  -  oo 

e«“‘’*rect(qa;„)  =  —  ^ 

?=0  °  k=-oo 


sin(f) 

% 


+  j 


OO 

/ 


sin(|)e 


‘-(ik-T) 


-  r) 


-dr 


(188) 


where  7*,  =  t  —  — . 

The  principal  value  integral  in  Equation  (188)  can  be  evaluated  using  residue  calculus.  Specif¬ 
ically,  the  integrand  is  first  broken  into  two  parts  using  Euler’s  identity  for  the  sine  function.  One 
then  allows  r  to  become  complex  where  each  integrand  part  has  different  regions  of  convergence. 
By  carefully  applying  Cauchy’s  theory,  one  can  relate  the  principal  value  integrals  to  the  sum  of  the 
half  residues  at  r  =  0  and  r  =  7*..  The  resulting  evaluation  gives 


^e^®“°*rect(ga)o)  = 

g=0 


OO 

“i  XI  i  >  0<a;o<l 


fc=— 00 


(189) 


-X 

Ulo  Z— / 


7k 


1  <  Wo 


fc=— 00 


which  can  be  verified  numerically.  Note  Equation  (189)  is  valid  for  all  values  of  t  and  positive  values 
of  Wo  as  long  as  Wo  7^  ^  for  integer  q.  This  restriction  is  a  consequence  of  the  “edge”  of  the  rect 
function  coinciding  with  one  of  the  sample  points  in  the  left  hand  side  of  Equation  (189). 

The  right  hand  side  of  Equation  (189)  provides  an  alternative  representation  which  becomes 
attractive  for  small  values  of  Wq.  In  that  case,  the  left  hand  side  will  require  many  terms  where 
the  corresponding  wide  separation  between  successive  values  of  ■jk  will  cause  the  right  hand  side 
summation  to  quickly  converge.  To  demonstrate  this  fact,  consider  the  particular  case  of  Wo  =  0.003 
andf  =  4.4429.  The  left  hand  side  summation  in  Equation  (189)  terminates  after  166  terms  resulting 
in  a  value  of  (60.299, 120.34).  Table  5  lists  the  partial  summations  of  the  right  hand  side  (RHS)  for 
an  increasing  number  of  terms.  The  relative  error  between  these  summations  and  the  left  hand  side 
(LHS)  is  also  given.  This  example  clearly  illustrates  the  inverse  bandwidth  nature  between  the  two 
representations  as  well  as  the  importance  of  such  a  transformation. 
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K 

K 

RHS 

k=-K 

0 

1 

2 

(60.198,120.47) 

(60.365,120.25) 

(60.282,120.36) 

0.00124 

0.00081 

0.00021 

Table  5.  Sequence  of  Partial  Sums  for  the  RHS 


138 


Appendix  B.  Summation  Acceleration  Algorithms 


When  using  the  Method  of  Moments  (MM)  to  solve  electromagnetic  scattering  problems, 
infinite  geometries  invariably  lead  to  infinite  summations.  The  semi-infinite  arrays  of  this  research 
are  no  exception.  In  seeking  a  numerical  solution,  one  is  forced  to  truncate  the  summations  in  one 
way  or  another.  The  obvious  choice  is  to  terminate  the  summation  after  a  given  number  of  terms,  but 
determining  how  many  terms  to  include  before  truncation  is  not  always  a  straight-forward  problem. 
As  is  often  the  case,  the  obvious  choice  is  not  always  the  “best”  choice.  This  appendix  describes 
several  algorithms  which  take  advantage  of  the  convergence  characteristics  of  the  infinite  summation 
to  accelerate  its  convergence.  It  turns  out  the  two  types  of  summations  used  in  this  research  (spatial 
and  spectral)  have  different  convergence  behaviors  which  affect  the  performance  of  each  algorithm. 
These  algorithms  are  the  e-algorithm,  the  circle-fit  algorithm,  and  the  p-algorithm. 

B.l  e-Algorithm 

The  e-algorithm  was  proposed  by  Wynn  [41]  as  an  efficient  alternative  to  Shank’s  transforma¬ 
tion  [30]  for  estimating  the  limit  of  a  sequence  of  partial  sums.  Unlike  Shank’s  transformation,  which 
requires  a  complete  recalculation  of  the  entire  sequence  for  each  estimate,  Wynn’s  e-algorithm  builds 
on  the  previously  calculated  estimates.  The  result  is  a  pyramid  of  numbers  whose  peak  contains  the 
best  estimate  to  the  infinite  sum’s  limit. 

The  entries  in  the  pyramid  are  obtained  from  the  so-called  em{‘)  function  operating  on  a 
sequence  of  partial  sums,  Sn-  The  pyramid  is  formed  by  implementing  the  following  equations 

^0  ( Sn )  —  Sn 

"  l«o(5.+i)  -  %(S.)1 

^m{Sn)  =  ^m-2{Sn+l)  +  J-  7^  \  /  rt  m 

where  e_i(S'„_|_i)  =  0.  Graphically  the  process  can  be  represented 


(190) 

(191) 

(192) 
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Note  that  odd  m  values  correspond  to  intermediate  results  which  are  necessary  for  the  next  estimate, 
^2m{Sn)  =  em{Sn),  where  6^(0  is  the  Shank’s  transformation.  Clearly  only  two  additional  partial 
sum  calculations  are  required  to  obtain  the  next  estimate.  The  complete  sequence  of  partial  sums 
does  not  have  to  be  recalculated  resulting  in  tremendous  computational  savings. 

As  in  all  asymptotic  series  evaluations,  care  must  be  taken  when  applying  them  to  a  particular 
summation.  Specifically,  the  e-algorithm  works  best  for  partial  sums  which  tend  to  oscillate  around 
their  convergence  value  (e.g. ,  the  spatial  summations  of  this  research).  When  applied  to  monotonically 
decreasing  series  such  as  those  generated  by  the  spectral  summations,  the  e-algorithm’s  performance 
degrades  somewhat.  Even  so,  the  author  found  it’s  acceleration  significant  enough  to  use  it  for  both 
the  spatial  and  spectral  domain  summations  encountered  in  this  research. 

B.2  Circle-Fit  Algorithm 

The  second  acceleration  algorithm  of  interest  also  applies  primarily  to  oscillatory  summations. 
The  circle-fit  algorithm  is  a  modification  of  the  spiral  average  method  proposed  by  Skinner  [33].  Both 
are  based  the  observation  that,  for  a  given  test  location,  the  consecutive  spatial  domain  partial  sums 
form  a  spiral  about  the  convergence  point  in  the  complex  plane.  This  phenomena  is  due  to  the  fact 
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that  the  distances  between  each  successive  array  element  and  the  test  point  asymptotically  approach 
the  inner-element  spacing.  By  estimating  the  center  of  the  spiral,  one  obtains  an  approximation  for 
the  eventual  convergence  point. 

The  primary  distinction  between  Skinner’s  spiral  average  algorithm  and  the  circle-fit  algorithm 
lies  in  the  technique  used  to  estimate  the  center  of  the  spiral.  Skinner  determines  the  number  of  partial 
sums  required  to  have  a  complete  phase  rotation  in  the  complex  plane.  He  then  performs  an  arithmetic 
average  of  the  these  partial  sums  to  obtain  his  estimate.  When  the  hmer-element  spacing  does  not 
lead  to  an  integer  number  of  phase  rotations,  Skinner  must  choose  an  error  parameter  and  search  for 
the  correct  number  of  partial  sums  required  to  obtain  a  phase  rotation  within  an  error  percentage  of 
an  integer  number  of  rotations. 

In  contrast,  the  circle-fit  algorithm  does  not  require  an  integer  number  of  phase  rotations.  The 
only  requirement  is  that  there  is  some  phase  rotation  over  the  collection  of  partial  sums.  As  the  name 
implies,  the  algorithm  simply  attempts  to  fit  a  circle  to  the  sequence  of  partial  sums.  The  center  of 
the  circle  provides  an  estimate  to  the  infinite  sum. 

As  an  illustration,  consider  the  two-dimensional  Green’s  function  for  a  semi-infinite  wire  array 

.  00 

G{x,y)  =lim  ^  V (193) 

g— >00  4  * 

5=0 

For  large  q  values  one  can  use  the  asymptotic  form  of  the  Hankel  function  to  express  the  summand 
as[l] 


s,  =  HP{py/{qd^  -  xf  4- 
~  Hf\l3\q\d^),  q-^oo 


2i 


TT, 


OO 


(194) 


Note  that  each  summand  contributes  a  component  with  a  constant  phase  increment  and  a  monotonically 
decreasing  amplitude.  Figure  117  depicts  a  typical  sequence  of  partial  sums. 
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Figure  117.  Circle-fit  Algorithm  for  Two-dimensional  Free-space  Periodic  Green’s  Function 

The  sequence  was  generated  using  Equation  (193)  with  the  parameters:  dx  =  0.21A,  x  =  0.0, 
y  =  0.01  A.  The  figure  also  contains  a  plot  of  the  “best  fit”  circle  to  the  last  five  partial  sums.  The 
center  of  this  circle  represents  the  current  estimate  of  the  infinite  sum.  To  obtain  a  final  estimate, 
one  compares  the  current  estimate  with  the  estimate  generated  from  the  next  five  partial  sums.  When 
the  difference  between  subsequent  estimates  reaches  a  specified  tolerance,  one  calls  the  process 
converged.  Note  in  this  case  the  estimates  are  generated  after  every  fifth  partial  sum.  In  general,  a 
new  estimate  could  be  generated  using  any  number  of  previous  partial  sums  (greater  than  three)  every 
time  a  new  partial  sum  is  calculated. 

The  method  for  fitting  the  circles  comes  from  a  paper  by  Moura  and  Kitney  [25].  In  it  they 
describe  how  to  determine  the  circle  which  minimizes  the  square  error  between  itself  and  the  set  of 
data  points.  Simply  stated,  the  problem  amounts  to  finding  a  center  point,  {xo,yo),  and  a  radius,  Vg, 
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such  that  the  error  defined  by 


^  K®"  ~  “  ^o] 


is  minimized.  This  is  accomplished  by  forcing  the  partial  derivatives,  ■§^,§^,  and  to  be  zero.  The 
latter  constraint  gives  a  relationship  between  the  data  set,  (a:*,  i/i),  the  circle’s  center  point,  {xg,  Vo), 


and  its  radius,  Vg 


From  the  remaining  constraints,  Moura  and  Kitney  derive  a  matrix  equation  whose  solution  gives  the 
circle’s  center,  [xg^yg). 

r Cl  1  r Gii  Ci2 ]  \ xg] 


G2I  C22  Vo 


(197) 


where 


N  ^  r  ^ 

i=l  L*=l 


N  N 


Gi2  (21)  -  4  ^  ^  Xj  ^  Vk 


y=i  k=i 


i=l  L*=l 

,  N  N 


N  N 


2  =  1  J  — 1  2=1  1;  =  1  2  =  1  2  =  1 


N  N 


N  N 


=  4e^/.’Ej'>-vE-?E!'‘+e-.“!'<+e*'.’ 


2  =  1  j  =  l 


2=1  fc=l  i  =  l 


Significantly,  this  method  is  a  direct  method^,  a  fact  crucial  in  maintaining  the  circle-fit  algorithm’s 
computational  efficiency. 


^  A  direct  method  implies  the  circle’s  defining  parameters  are  calculated  from  the  data  points  in  one  step.  There  is  no 
iterative  process  to  refine  an  initial  guess. 
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Relative  Error 

Raw  Sum  (#  terms/cpu  secs) 

Circle-fit  Algorithm 

e-Algorithm 

l.OE-1 

10/0.001 

5  /  0.001 

2  /  0.001 

l.OE-2 

977  /  0.4 

15/0.004 

4/0.002 

l.OE-3 

24834/50.7 

45/0.011 

6  /  0.003 

l.OE-4 

NC 

190  /  0.042 

6  /  0.003 

l.OE-5 

NC 

880/0.167 

8  /  0.004 

l.OE-6 

NC 

3865  /  0.766 

10  /  0.005 

l.OE-7 

NC 

11380/2.375 

12/0.005 

Table  6.  Circle-fit  and  e  Algorithms  Applied  to  Equation  (193) 


Table  6  clearly  demonstrates  this  efficiency  by  comparing  the  convergence  characteristics  of 
the  two  acceleration  algorithms  to  the  unaccelerated  sum  of  Equation  (193).  The  parameters  here  are 
the  same  as  those  used  to  generate  Figure  1 17’s  plot.  In  the  table  the  letters  “NC”  indicate  the  sum  did 
not  converge.  Convergence  is  determined  by  calculating  the  relative  error  between  each  algorithm’s 
estimate  and  a  “truth”  value  obtained  by  setting  the  convergence  tolerance  to  machine  accuracy  for 
either  algorithm.  The  relative  error  is  then  defined  by 


relative  error  = 


estimate  -  truth 
truth 


(198) 


Also,  the  CPU  times  are  those  obtained  on  a  Sun  SparclO  workstation.  Based  on  these  results,  the 
author  chose  the  e-algorithm  over  the  circle-fit  algorithm. 


B.3  p-Algorithm 

The  final  algorithm  considered  reportedly  performs  best  on  monotonically  decreasing  sum¬ 
mations.  The  p-algorithm  was  proposed  by  Singh  and  Singh  as  a  more  efficient  alternative  to  the 
e-algorithm  [31].  Originally  derived  by  Wynn  [42],  the  algorithm  can  be  computed  by  the  following 
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Relative  Error 

Raw  Sum  (#  terms/cpu  secs) 

e-Algorithm 

p-Algorithm 

l.OE-1 

5  /  0.003 

4/0.003 

2/0.003 

l.OE-2 

60/0.015 

14/0.005 

2  /  0.003 

l.OE-3 

607/0.129 

132/0.056 

4/0.003 

l.OE-4 

6078/  1.4 

136/0.056 

4/0.003 

l.OE-5 

NC 

730/0.123 

6/0.003 

l.OE-6 

NC 

NC 

6/0.003 

l.OE-7 

NC 

NC 

8  /  0.003 

Table  7.  e  and  p  Algorithms  Applied  to  Equation  (202) 


equations 


PO  (  )  — 

Pl{Sn)  = 

Pm{Sn)  — 

where  p-i{Sn+i)  =  0.  Notice  the 
of  the  last  equation.  It  is  claimed 
over  the  e-algorithm  when  operating  on  a  monotonically  decreasing  summand. 

As  a  demonstration,  consider  the  infinite  summation 

OO  - 

s  =  E  p  (202) 

Jb=l 

Table  7  contains  a  comparison  of  the  e-algorithm  and  p-algorithm  along  with  the  unaccelerated  sum. 
Based  on  this  data,  the  p-algorithm  would  be  the  algorithm  of  choice.  In  fact  in  their  paper,  Singh 
and  Singh  present  relative  error  plots  of  this  algorithm  applied  to  both  spectral  and  spatial  Green’s 
functions  similar  to  those  used  in  this  research  which  clearly  indicate  a  performance  advantage  over 
the  e-algorithm. 

Unfortunately,  the  author  was  not  able  to  obtain  the  same  performance.  Table  8  shows  the 
algorithm  applied  to  the  discrete  spectral  sum  of  Equation  (203)  with  parameters:  =  0.21A, 


[po(‘5^n-l-l)  Po('^n)] 

Pm—2{Sn+l)  “b 


m 


[Pm— l('S'n+l)  Pm— l{Sn)] 


(199) 

(200) 

(201) 


similarity  to  the  e-algorithm.  The  difference  lies  in  the  numerator 
this  small  change  gives  the  p-algorithm  a  tremendous  advantage 
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Relative  Error 

Raw  Sum  (#  terms/cpu  secs) 

e-Algorithm 

p-Algorithm 

l.OE-1 

21  /  0.009 

8  /  0.004 

4/0.004 

l.OE-2 

70  /  0.024 

14/0.006 

16/0.005 

l.OE-3 

130/0.042 

20/0.008 

136/0.075 

l.OE-4 

196  /  0.064 

32/0.012 

170/0.109 

l.OE-5 

264/0.087 

52/0.018 

392  /  0.473 

l.OE-6 

333/0.106 

74/0.028 

NC 

l.OE-7 

404/0.127 

84/0.036 

NC 

Table  8.  e  and  p  Algorithms  Applied  to  Equation  (203) 


X  =  0.0,  y  =  O.OOIA. 


=  to  E 


(203) 


The  p-algorithm  actually  requires  more  terms  than  the  raw  sum.  Based  on  these  results  and  other 


numerical  experiments,  the  author  decided  the  e-algorithm  provided  the  most  reliable  acceleration  in 
both  the  spatial  and  spectral  domains.  This  fact  coupled  with  the  algorithm’s  relative  simplicity  made 
the  e-algorithm  the  method  of  choice. 
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Appendix  C.  Continuous  Spectrum  Integral:  Uniform  Asymptotic  Expansion 

The  application  of  the  Method  of  Moments  (MM)  to  infinite  periodic  structures  results  in 
infinite  summations.  The  summation  acceleration  algorithms  of  Appendix  B  provide  one  approach 
for  evaluating  these  summations.  Another  approach  involves  transforming  the  summand  to  another, 
more  convergent,  form.  The  one-sided  Poisson  sum  used  in  this  research  is  an  example  of  this 
technique.  Unfortunately,  the  one-sided  Poisson  sum  generates  a  continuous  spectrum  integral  which 
is  often  difficult  to  evaluate  numerically.  In  an  attempt  to  mitigate  this  problem,  the  author  investigated 
the  following  asymptotic  integral  evaluation.  Although  not  universally  applicable,  the  asymptotic 
integral  evaluation  significantly  shortens  the  computation  time  required  for  many  of  the  integrals 
encountered  in  this  research. 

The  continuous  spectrum^  integral  has  the  general  form 

Ik  =  dv  (204) 

c 

where  (7  is  an  infinite  contour  in  the  complex  p  plane,  q{p)  is  an  analytic  function,  and  f{p)  has  a 
simple  pole^  located  on  C.  The  parameter  Q,  is  positive  and,  most  important,  often  large.  When  this 
occurs,  one  can  use  an  asymptotic  approximation  for  and  still  maintain  reasonable  accuracy. 

At  its  core,  the  asymptotic  integral  evaluation  is  based  on  the  exponential  behavior  of  the 
integrand.  Suppose  Re{q{p)}  in  Equation  (204)  has  a  local  maximum  Pg,  known  as  a  saddle  point, 
somewhere  along  C.  In  the  neighborhood  of  Pg,  Re{q{Ps)}  >  Re{q{p)].  Since  is  a  large  positive 
value,  the  magnitude  of  the  exponential  function,  |  | ,  will  be  greatest  at  Pg  and  decay  very  rapidly 

away  from  Pg.  As  long  as  f{p)  is  analytic  in  the  neighborhood  of  Pg,  one  can  approximate  R  with 
only  the  integral  contribution  from  the  local  maximum. 

^This  form  is  actually  the  angular  spectrum  resulting  from  the  transformation  px  =  sin(t/) . 

^For  the  non-free- space  geometries,  /(v)  may  also  contain  singularities  associated  with  surface  waves  trapped  in  the 
stratified  media. 
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At  this  point,  one  usually  applies  a  little  finesse  in  an  effort  to  make  the  local  evaluation  as  easy 
as  possible.  For  example,  if  /(z/)  is  slowly  varying  in  the  neighborhood  of  Ug,  it  can  be  approximated 
a)  and  pulled  out  of  the  integral 

4~/K)  J  n->oo  (205) 

\v-v,\<f 

One  can  then  perform  a  Taylor  expansion  of  q{u)  about  Ug.  Retaining  the  first  two  non-zero^  terms, 
one  arrives  at  the  “canonical”  form 

h J  Q^oo  (206) 

|l/— I/a|<e 

Since  the  integrand  is  exponentially  dominant,  one  can  extend  the  range  of  integration.  Then,  with 
the  appropriate  variable  substitution,  one  can  express  the  integral  in  terms  of  the  gamma  function. 
Closed  form  evaluation  gives  the  final  asymptotic  form 

Unfortunately,  the  integrals  encountered  in  the  dissertation  rarely  display  the  desired  behavior 
without  some  manipulation.  For  example,  one  must  deform  C  into  what  is  commonly  called  the 
steepest  descent  path,  Csdp-  Along  Csdp,  the  integrand  displays  the  desired  sub-dominant  behavior 
away  fi-om  the  saddle  point.  Another  difficulty  occurs  when  /(z^)’s  singular  point  lies  in  the  vicinity 
of  the  saddle  point.  One  must  then  decompose  f[v)  into  a  regular  and  singular  part.  The  regular  part 
leads  to  the  gamma  function  as  before.  The  singular  part  leads  to  error  function  or  Fresnel  integrals 
which  are  well  tabulated.  Felsen  and  Marcuvitz  provide  an  excellent  overview  of  this  process  in  their 
chapter  on  the  asymptotic  evaluation  of  integrals  [11]. 

^Since  Vg  is  a  maximum,  q'{vg)  =  0  where  the  prime  indicates  differentiation  with  respect  to  the  argument.  Thus  the 
next  non-zero  term  is  q"{va). 
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To  illustrate  the  benefits  of  asymptotic  evaluation,  consider  the  integral  arising  from  the  two- 
dimensional  wire  array  in  free-space 


h 

I 


sin(z/)  -  pj, 


dv 


(208) 


where 

k\ 

=  +  ^  (209) 

dx 

and  C  is  shown  in  Figure  6.  The  parameter  ppt  represents  the  electrical  distance  from  a  point  one 
half  the  inner-element  spacing  to  the  left  of  the  array  edge  and  the  angle  measured  from  the  array 
normal.  Thus  for  many  of  the  terms  in  the  coupling  matrix,  /3pt  is  a  “large”^  parameter. 

To  perform  an  asymptotic  evaluation,  one  must  distort  the  principle  value  contour  to  the  steepest 
descent  path,  accounting  for  any  pole  encirclement  along  the  way.  In  addition,  the  principle  value 
contour  must  be  completed  in  order  to  invoke  Cauchy  residue  theory  to  distort  the  contour.  One  can 
incorporate  both  factors  in  the  analysis  through  the  combined  residue  contribution 


Ik  ICsDP  I  residue 


(210) 


Figure  118  depicts  these  contours  for  a  typical  saddle  point/grating  lobe  pole  configuration. 
The  grating  lobe  pole  occurs  at  Vo  =  sin“^(pa;)  and  the  saddle  point  at  Vg  =  The  steepest  decent 
path  is  defined  by  Im{q{v)}  =  Im{q[vs)'\  which  intersects  the  real  axis  at  Vg  and  one  of  the  vertical 
legs  at  i/i  =  ±  cosh“^  (|8in((^0|)  ’  <i®P®i^ding  on  the  sign  of  Vg. 

Felsen  and  Marcuvitz  provide  a  uniform  asymptotic  evaluation  valid  even  when  the  grating 
lobe  pole  and  the  saddle  point  come  arbitrarily  close  to  each  other.  Expressed  in  terms  of  the 

■^Of  course,  the  definition  of  “large”  depends  on  the  accuracy  one  desires. 
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Figure  118.  Angular  Spectral  Integral  Steepest  Descent  Path 
complementary  error  function,  the  evaluation  is 


SDP 


~  e 


I  27r 

I3pt 


±  j2  -e 


_ p -iPpi  [cOs( Vo-^t) 


|^sm(0t)  p^y^2[cos(^'o  -  0*)  -  1]J 

Q  ^ l3pt[co^{uo  -  ^t)  -  1]^  I 


(211) 


where 


OO 


Q{x)  ~  j  ^ 

(212) 

X 

Py  =  Vl 

(213) 

Note  that  py  is  defined  to  be  either  positive  real  or  negative  imaginary.  Abramowitz  and  Stegun  [1] 
give  efficient  methods  for  evaluating  Q{x)  with  a  complex  argument  without  resorting  to  numerical 


integration.  One  can  then  calculate  the  remaining  half-pole  contribution  from 


^residue  —  i 


jTte 


-iPptCOs(Vo-<j>t) 


Py 


(214) 


To  obtain  the  correct  residue  contributions  one  must  choose  the  proper  signs  in  Equations  (211)  and 
(214).  Specifically,  the  upper  signs  are  chosen  if  one  of  the  following  conditions  are  satisfied. 


1.  Vg  >  0  and  Re{i/o}  <  Vg 

2.  Vg>  ^  and  Im{vo]  >  Vi,  or 

3.  Vg  <  0,  Re{vo]  <  Vg,  and  Im{vo]  >  —Vi 
If  none  are  true,  then  one  chooses  the  lower  signs. 

As  in  any  approximation,  the  question  of  accuracy  invariably  arises.  Figure  119  attempts 
to  provide  an  answer  to  this  question  by  comparing  the  asymptotic  and  numerical  evaluations  of 
Equation  (208)  using  the  parameters  ^  s*  =  —0.866,  =  O.IA  and  k  =  0.  The  figure  plots 

the  relative  error  magnitude  defined  by 


relative  error  = 


asymptotic  —  numerical 
numerical 


(215) 


as  a  function  of  0*  for  various  values  of  pt.  Since  the  test  locations  in  this  dissertation  are  always 
a  wire  radius  out  of  the  array  plane,  it  is  clear  stays  close  to  -1-90°  for  most  integral  evaluations. 
Surprisingly,  pt  can  be  as  small  as  O.IA  over  a  wide  ^t  range  while  still  maintaining  a  relative  error 
less  than  10%.  This,  coupled  with  the  computational  savings®  shown  in  Table  9,  make  the  asymptotic 
evaluation  a  very  attractive  alternative  to  a  numerical  solution. 


®The  CPU  times  are  averages  over  91  integral  evaluations.  The  evaluations  were  performed  on  a  Sun  SparclO 
workstation  using  the  MATLAB  software  package. 
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Relative  Error  Magnitude 


Figure  119.  Numerical  versus  Asymptotic  Integral  Evaluation  Error 


P(A) 

Asymptotic  (cpu  secs) 

Numerical  (cpu  secs) 

0.1 

0.06 

9.08 

0.5 

0.06 

9.38 

1.0 

0.05 

9.48 

Table  9.  Asymptotic  and  Numerical  Integral  Evaluation  Time  Comparison 
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Appendix  D.  Array  Scanning  Method 


Consider  the  problem  of  a  single  radiating  element  in  the  presence  of  a  stratified  dielectric 
medium.  One  must  be  able  to  determine  the  fields  from  such  an  element  in  order  to  calculate  the 
coupling  matrices  found  in  this  research^ ,  Typically  one  might  use  a  non-free-space  Green’s  function 
to  account  for  the  slab  boundaries.  Unfortunately,  this  type  of  Green’s  function  often  produces 
difficult  to  evaluate,  highly  oscillatory  integrands.  This  appendix  briefly  develops  an  alternative 
method  called  the  Array  Scanotiing  Method  (ASM).  The  author  uses  a  two-dimensional  wire  problem 
to  illustrate  the  ASM,  thereby  avoiding  a  tremendous  amount  of  mathematical  baggage  which  can 
mask  the  method’s  simplicity. 

Figure  120  depicts  a  single  infinite  wire  buried  in  a  dielectric  slab.  One  begins  by  expressing 
^The  coupling  matrix  entry  is  proportional  to  the  field  radiated  by  the  basis  function  at  tbe  test  location. 


the  field  due  to  the  wire  in  a  homogeneous  media  as 


(216) 


where  T  =  Since  one  has  total  freedom  in  defining  the  array,  one  can  choose  the 

inner-element  spacing  to  confine  the  integration  range  to  the  interval  [—1,1].  This  corresponds  to 
scanning  the  array  phasing  through  an  angle  range  restricted  to  real  space,  hence  the  name  “array 
scanning”  method. 

At  this  point,  one  can  express  the  field  from  the  array  in  the  spectral  domain  via  the  Poisson 
sum  formula.  It  remains  only  to  introduce  the  T-factor  concept[26,  28]  to  incorporate  the  dielectric 
slab  into  the  problem.  The  T-factor  will  not  alter  the  phase  relationship  between  the  array  elements 
and  therefore  does  not  alter  the  scan  integration.  Thus,  one  is  left  with  an  expression  for  the  field 
from  a  single  wire  in  terms  of  the  field  from  an  infinite  array  of  wires. 
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One  may  wonder  what  has  been  gained  through  this  series  of  manipulations.  It  appears  one 
has  exchanged  a  straight-forward  integral  evaluation  for  an  infinite  sum  of  integral  evaluations.  The 
ASM  derives  its  efficiency  from  the  exponential  damping  characteristics  of  the  evanescent  spectral 
terms.  Thus  only  a  few,  much  simpler,  integral  evaluations  are  required  for  convergence.  The  only 
penalty  occurs  when  the  T-factor  contains  surface  wave  singularities.  In  this  case,  the  most  expedient 
remedy  involves  introducing  a  slight  loss  in  the  dielectric  and  then  integrating  finely  near  the  former 
singular  point. 
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Appendix  E.  Physical  Basis  Function  (PBF)  in  the  Excitation  Vector 

Throughout  this  dissertation,  the  author  treats  the  PBF  coefficient  as  an  unknown  quantity. 
Recall  the  PBF  was  introduced  as  an  <3  priori  encapsulation  of  the  “infinite”  nature  of  the  central 
portion  of  a  large  array.  One  might  ask  why  not  go  one  step  further  and  set  the  PBF  coefficient  value 
to  that  determined  by  a  full  infinite  array  analysis.  It  could  then  be  included  in  the  excitation  vector, 
reducing  the  matrix  dimension  and  eliminating  the  question  of  how  to  test  the  PBF.  This  appendix 
documents  the  author’s  investigation  of  this  idea. 

Consider  the  simple  problem  of  a  large  wire  array  in  a  homogeneous  media  shown  in  Figure  121. 
Breaking  the  problem  into  two  semi-infinite  arrays,  one  could  formulate  the  subproblem  Moment 
Method  (MM)  matrix  equation  in  the  following  two  ways.  Treating  Jq  as  an  unknown,  one  has  from 


O 


"Finite"  Portion  "Infinite"  Portion  "Finite"  Portion 

(100  Elements)  (400  Elements)  (100  Elements) 


<t)‘=  200  ° 
d  =  0.3  X 


Figure  121.  600  Element  Wire  Array  in  Homogeneous  Media 
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Equation  (16) 


Blip,)  =  - 


UjJ, 

T 


B 


Jo  Y,  -  p,\)e-’>’‘  l‘<  +  Y,  ■’l■Ho‘\P\P^  ~  ft!) 


g=0 


b=l 


(219) 


where  the  right  hand  side  contains  the  unknown  current  coefficients.  Recall  'p^  is  the  test  function 
location,  pj  is  the  basis  function  location,  and  is  the  PBF  element  location.  The  second  formulation 
treats  the  PBF  coefficient,  Jq.  as  a  known  resulting  in 


Bi(P,)  +  ^  =  -X  E  -  PJ)  (220) 


B 


q=0 


6=1 


One  can  determine  Jq  from  a  full  infinite  array  formed  about  the  PBF  test  location,  pp. 


Jo  =  - 


iBUPo) 


hPWPo  -  P,\)e 

q=—oc> 


(221) 


The  following  four  figures  compare  the  currents  calculated  with  these  two  formulations  (denoted 
by  x’s)  to  those  calculated  from  a  600  unknown  MM  analysis  (denoted  by  o’s).  Specifically,  the 
current  amplitudes  for  the  first  110  edge  elements  are  shown.  Figures  122  and  123  plot  the  currents 
determined  by  Equation  (219)  with  100  unknown  edge  elements.  For  clarity,  the  first  ten  PBF 
elements  are  also  shown.  Note  how  the  edge  elements  near  the  PBF  are  clamped  to  the  PBF  value. 
Now  consider  Figures  124  and  125  where  Equations  (220)  and  (221)  are  used.  In  this  formulation, 
the  edge  perturbations  near  the  PBF  are  no  longer  mitigated. 

To  understand  why  this  is  so,  consider  the  matrix  equation  representing  Equation  (219). 


Zoo 

Zob 

Jo 

.  . 

Zto 

Ztb 

.  . 

(222) 
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where  the  0  subscript  indicates  the  PBF.  Compare  Equation  (222)  with  that  of  the  second  formulation. 


K  -  4Zto 


(223) 


where  one  calculates  Jg  from  a  fully  infinite  array.  Even  if  Jg  =  Jg,  Equation  (223)  does  not  contain 
the  constraint 

Fg  =  JgZgg  +  [Zgfc]  [J']  (224) 

Thus  the  second  formulation  is  missing  some  of  the  coupling  information  between  the  PBF  and  the 
edge  elements. 

Based  on  the  previous  figures,  it  is  apparent  this  coupling  information  plays  a  critical  role  in 
correcting  the  edge  element  currents  nearest  the  PBF.  As  a  result  of  this  investigation,  the  author 
determined  it  is  better  to  treat  the  PBF  coefficient  as  an  unknown  despite  what  may  appear  to  be  an 
unnecessary  increase  in  the  resulting  matrix  size.  One  can  always  compare  the  PBF  coefficient  to  that 
obtained  from  an  infinite  analysis  “after-the-fact”  to  check  if  a  sufficient  number  of  edge  elements 
were  included  in  the  calculation. 
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Figure  122.  Left  1 10  Elements  Based  on  First  Formulation  (Leading  Edge  of  Array  is  Element  1) 


Figure  123.  Right  110  Elements  Based  on  First  Formulation  (Trailing  Edge  of  Array  is  Element 
600) 
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Figure  124.  Left  1 10  Elements  Based  on  Second  Formulation  (Leading  Edge  of  Array  is  Element  1) 


Figure  125.  Right  110  Elements  Based  on  Second  Formulation  (Trailing  Edge  of  Array  is  Element 
600) 
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